A QUANTITATIVE DIFFERENTIAL EQUATION 
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MALWINA LUCZAK 

Abstract. We consider a Markovian load balancing model on a fully- 
connected network, where calls have Poisson arrivals and exponential 
durations. The endpoints of each call are uniform over all the links of 
the network. Each call is routed either along the link connecting its 
endpoints, or, if the direct route is unavailable, along a two-link path 
between them, via an intermediate node. We use an explicit and simple 
coupling to show a strong concentration of measure property, and deduce 
that the evolution of the process may be approximated by a differential 
equation. The technique is likely to be useful in other settings. 



S : intro 



1. Introduction 

We consider a class of online routing problems in continuous time, where 
calls have Poisson arrivals and exponential durations, studied earlier in [l|; 

MMMM- 

The setting is as follows. For each n € N, we have a fully connected 
communication graph Kn, with node set Vn = {1, ■ ■ ■ ,n} and link set En = 
{{u,v} : 1 < u < V < n}. Each link {u,v} £ En has capacity C = 
C{n) < oo, where C{n) G Z+ for each n. Each arriving call is to be routed 
either along a link {u, v} or along a path between u and v consisting of a 
pair of links {u,w} and {v,w}, for a pair u,v of distinct nodes (endpoints 
of the call) and some intermediate node w ^ u,v, if possible. A call in 
progress will use one unit of capacity of each of the links it occupies, for 
its entire duration. Calls arrive as a Poisson process with rate A (2), where 
A = A(n) > 0. The endpoints of each call are uniform over the links of 
the complete graph Kn- If a call is for nodes u and v, then we route it 
over the direct link {n, v} between u and v if possible, that is if {u, v} 
has fewer than C calls currently using it. Otherwise, we pick an ordered 
list of d = d{n) possible intermediate nodes {wi, . . . ,Wci) from Vn \ {u,v}, 
uniformly at random with replacement, and the call is routed along one 
of the two-link routes {{u,wi},{v,wi}}, . . . ,{{u,Wii}-,{v,Wii\}-, chosen to 
minimise the larger of the current loads on its two links, subject to the 
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capacity constraints. Ties are broken in favour of the first 'best' route in the 
ordered hst. If none of the d two-hnk paths is available, then the call is lost. 
Call durations are unit mean exponential random variables, independent of 
one another and of the arrivals and choices processes. 

Here, we focus on the analysis of this algorithm as n tends to infinity. 
We prove that, asymptotically, for suitable initial conditions and suitable 
functions A(n), d(n) and C{n), for each node v, the proportion of links 
at V that carry k calls is well approximated by the solution of a differential 
equation. (Note that, when A, d and C vary with n, there is no single limiting 
differential equation, but rather a sequence of approximating differential 
equations, with dimension tending to infinity if C ^ oo with n.) 

Such law of large numbers results have been difficult to prove in this and 
related models, due to apparent and potentially strong dependencies be- 
tween system elements (in this context, links). Here we are able to prove 
that these dependencies are negligible; it turns out that, in a suitable sense, 
links in certain collections evolve approximately independently of one an- 
other. Our technique appears to be new, and is likely to be useful in other 
settings. It relies on a coupling, which is used to prove that slowly-changing 
functions of the process (for instance, the number of links around a node v 
with load exactly k, for each node and each k G {0, 1, . . . , C}) are well con- 
centrated at each time t. Thanks to the strong concentration of measure, it 
is then possible to show that the expected drifts of functions of interest fac- 
torise approximately, leading to a differential equation approximating these 
functions. The basic principle of our approach is, in essence, rather simple; 
however, there are considerable complications arising from the complicated 
evolution of the process in question. 

For each link e = {u, v} E E^, let X^"^ (e, 0) denote the number of calls in 
progress at time t which arc routed along the link e, that is the number of 
directly routed calls between the end nodes u and v of e that are in progress 
at time t. For each link e = {u, v} G En and node u) G \ e, let X^"'\e, w) 
denote the number of calls in progress at time t which are routed along the 
path consisting of links {u,w},{v,w}, that is the number of calls between 
the end nodes u and w of e routed via w that are in progress at time t. We 
can Xj^"^ = {X^'^\e,0),X^'^\e,w) : e e En,w e Vn \ e) the load vector at 
time t, and let S = {0, 1, . . . , C}"^"^^) denote the state space, containing 

the set of all possible load vectors. Then = (Xj"'^)t>o is a continuous- 
time discrete-space Markov chain. We will normally drop the superscript n, 
to avoid unnecessarily cluttering the notation. 
Let S' C be the set of load vectors x such that 



x\\i = 




that is the subset of the state space consisting of those load vectors where 
the total number of calls in the network is at most 6 A (2). Let T be the 
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first time t such that Xt S. Then T is a stopping time with respect to 
the natural filtration of process X. Let X'^ be the Markov process given 
as follows: Xf = Xt for t < T, and Xf = Xt for t>T. In what follows, 
we omit the superscript T, again to lighten notation, and write simply X, 
but it is to be understood that we only consider the evolution of the process 
until it leaves S. In the case where the link capacity C is constant, it would 
be slightly simpler to proceed without restricting the evolution of X to S, 
but we will not bother treating this case differently. 

Given a load vector x & S and a link e = {u, v} G E^, let x{e) denote the 
load of link e. Then 

x({u,v}) = x{{u,v},0) + {x{{u,w},v) + x{{v,w},u). 

Given a load vector x, node v and k G Z"^, let fv,k{x) be the number of links 
{v,w}, w V, in X such that x{{v,w}) = k (that is, the number of links 
with one end v carrying exactly k calls). 
For a vector ^ = (^(A;) : = 0, . . . , C), let 

k=0 

Define F : R'^+^ R'^+^ by 

FkiO = - I) - >^ak) + Xgk-iiO - >^9ki0 - k^k) 

+ + + 1), 0<A;<C, 

FciO = XaC-l) + X9C-i{0-mC), (1-1) 
where functions gj, j = 0, . . . ,C — 1, are given by 

9j{0 = nc)cm{< - a< jfY-\i - a< J - If)'-' 

r=l 

i=j+l r=l 

(i-e(<j-i)')"-(i-e(<jr)" 



i=j+i 



e(<jT-e(<j-i)2 
{<i-ifY-{i 



When d=l, 5,(0 = 2e(C)(l-C(C))e(j) for < j < C- 1. When d = 2, 
for < i < C - 1, 

= 24(C)C(J)C(<J)(2-(C(<J))'-(^(<J-1))') 
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j<i<C 



m i))'-(e(<^-i))'). 



We will show (see Lemma l7. ID that, for each d > 1, the function -F(^) is 
Lipschitz, with constant 8(P{X + 1)(C + 1)^, over the set of non-negative ^ 

such that X^fc=o ^(^) — ^' '^ith respect to the £00 norm. Hence we will prove 
that the equation 



has a unique solution starting from ^0 such that ^o(j) ^ for all j and 
Z^^o ^oO) = 1) valid for all times and such that, for all t, £,t{j) > for each 

j and EjLo^tO') = 1- 

In equation (jl.3p . the linear terms account for the effects of calls routed 
directly, and for call departures. Each function gj is proportional to the rate 
of arrivals of alternatively routed calls onto links that carry j calls. 

Given a pair of nodes u,v and j £ {0, . . . ,C}, let liv : 5 — )■ {0, 1} be 
defined by liv{x) = 1 if x{{u,v}) = j and iLvix) = otherwise. Thus liv 
is the indicator of the set of load vectors where the load of link {u,v} is j. 
Note that lu^ = liu, and that 1^,, is identically for each v and j. 
Let Si be the set of all states x such that ||x||i < '^XiJ^)- 

: S -^Rhe defined by 



Let functions 



max max 

u,v:u^v j,k 
1 



n 



(n-2)2 4^ 



X) E 1 



k 

vw' 



{x] 



(L4) 



1 



max max — 



:\fu,jix) - fv,j{x)\ 



max max — 

u,v:uj^v j n ■ 



/ J '^vw 



max max — 



4>\x) 



1 



max — 



E x{{u,v},w). 



(L5) 
(1-6) 



Let (/) = max{i;^^, cf^jcf)^}. This function is related to the function called 
by Crametz and Hunt in 01- 

We shall prove that, if (j){^o) is small, for instance (/>(Xo) = O^^^^ 

then (p{Xt) remains small for a time interval of order 1, and over that period 
each function (n — fvj(Xt) is well-approximated by the solution to the 
differential equation (jl.3p with initial condition ^o(j) = C'^ — l)~^/t) j(-^o) for 
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thm . main-result 



thm . main-result-d=l 



j = 0, . . . ,C. Note that, if </>^(Xo) is small, then, for each j, all the functions 

/„j(Xo) for different v £ Vn are nearly equal, so all the (n — 1)"-^/^, 

can be approximated by the same solution to the differential equation (jl.3p . 



Theorem 1.1. Suppose that A and are positive reals, and d and C posi 
tive integers. Set 7 = l/{2'^'^d^{8Xto + i)3e800dAto)^ ^nd suppose that 

^ ^,10\X + imC + lftle'/^ 



"^"^""VA'SAto 

Let Co satisfy ^o(i) > for all j = 0,..., C, and Yl^o Co(i) = 1- Let {^t) 
be the unique solution to the differential equation on [0,to]; subject to 

initial condition ^o- Let Xq be in Si. Let be the event that, for each 
V S Vn, each k £ {0, . . . ,C} and each t E [0, to]; 

\f,,k{Xt) -{n- < (sup fu,j{Xo) - (n - l)Co(i) 



+ 46(A + l)(to + l)d\C + 1)=^ (ncPiXo) + Sn^/^ 



; n 



^216{X+l)cfi{C+iyHo 



Then F{An) < e-hi°s'". 

In particular, suppose that there is vq £ Vn such that $,o(j) = :^;^fvo,j{^o) 
for j = 0, ...,C. Let A'^ be the event that, for each t < to, each k £ 
{0, . . . , C}, and each v G Vn, 

\f.,k{Xt) - (n - l)it{k)\ < 46(A + l)(to + l)d\C + l)^ 



X f 2n0(Xo) + 3n^/2 log n) e 



„216(A+l)d2(c+l)3to 



ThenF{A'J < e'i^^"^ 

For the special case d = 1, we obtain sharper bounds, replacing the term 
(C + 1)3 in the exponent by (C + 1). 

Theorem 1.2. Suppose that A and to are positive reals, and C a positive 
integer, and suppose that d = 1. Set j = l/{2^'^{8Xto + l)3e^°°^*o), and 
suppose that 



n > max — 



1 1 



. A 8Ato 

Let Co satisfy CoU) > for all j 



,10^(A + 1)4(C + 1)6^2, e8/7 
c 

j=0 



0,...,C, andZtoUj) 



1. Let (6) 

be the unique solution to the differential equation il.3\) on [0, to], subject to 
initial condition Co- Let Xq be in Si. Let Bn be the event that, for each 
V £ Vn, each k £ {0, . . . ,C} and each t £ [0, to], 



\f.,kiXt) -in- l)Ctik)\ < sup fuAXo) -in- l)Coij) 



+ 46( A + 1) (to + 1) (C + 1)3 (^n,/.(Xo) + Sn^/^ log 



; n e 



216(A+l)(C+l)to 



Then P(5„) < e 
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Suppose that, for each n, Xq ' = Xq a.s. for some deterministic load 
vector Xq"'* such that, for some constant c, (f>{x^^) < ^1^" and max^j \ {n — 

^)^^ fv,j{xl^^) — Co{j)\ ^ '^1^" ■ Suppose also that, as n — oo, X and to are 
bounded away from 0, and that XcfC'^to = o(logn) and dXto = o(loglogn). 
Then, for sufficiently large n, the condition on n in Theorem II. II is satisfied, 
and the theorem implies that, for e > 0, if is the event that, for each 
V £ Vn, each k £ {0, . . . , C}, and each t £ [0, to], 

|/,,fc(Xt)-(n-l)6(A;)|<nV2+^, 

then P(^^) — )• as n — )• oo. For d = 1, the corresponding conditions are 
that, as n — )• oo, A and to are bounded away from 0, and that XCIq = o(log n) 
and XIq = o(loglogn). 

In the simplest case where d = 1, if the direct link is at full capacity, 
only one two-link alternative route is considered, and it is used if there 
is spare capacity on both links. This case, with constant arrival rate A 
and constant capacity C, was first studied by Gibbens, Hunt and Kelly [1] 
and then by Crametz and Hunt [3], and Graham and Meleard 0]. For 

= 0, 1, . . . , C, let Y^'^\k) denote the proportion of links that carry k calls 
at time t in a system with n nodes. It is conjectured in and shown 
in [3] that, under suitable conditions, {Y^'^\k) : k = 0, ...,(7) converges 
in distribution as n — )• oo to a deterministic vector {Yt{k) : k = 0, . . . , C) 
obtained as the solution to the differential equation derived from the average 
drift of {Y^'^\k) : k = 0, ...,C), with appropriate initial conditions. The 
convergence result in is non-quantitative. Graham and Meleard do 
give a quantitative result concerning independence of small collections of 
links under the assumptions that initial link loads are iid and that initially 
there are no alternatively routed calls in the network. This result can be 
used to deduce a quantitative law of large numbers. (Also, it would be 
possible to quantify convergence in the more general case they consider.) 

In the case of A and C constant, and d = 1, Theorem ll.2l is a more refined, 
quantitative, version of the law of large numbers in [3|. Also, our result in 
this case is related to those in . Unlike [1] , we do not need to assume that 
initially all the nodes are exactly exchangeable. Instead, our law of large 
numbers result holds for a large class of deterministic initial states, and 
holds simultaneously for all nodes. Theorem 1 1 . 1 1 and the remaining cases in 
Theorem 11.21 are completely new. 

For d > 2 constant and constant A, this model is a variant of one that 
has attracted earlier interest. Luczak and Upfal study a version (both 
in discrete and in continuous time) where the total capacity of each link 
{u, v} is divided into three parts, one for 'direct' calls, one for indirectly 
routed calls with one end u and one for indirectly routed calls with one end 
V. Equivalently, each 'undirected' link {u,v} has capacity Ci(n) and is a 
first-choice path for calls between u and v; also, for each link {u, v} there 
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are two directed links, uv and vu, each with capacity C2{n). Link uv is used 
for indirectly routed calls with one end u and link vu is used for indirectly 
routed calls with one end v. The results of for the discrete-time model 
were strengthened and extended by Luczak, McDiarmid and Upfal [s*], who 
also studied the discrete-time version of the model that is the focus of this 
paper. The long-term behaviour of the continuous-time model was analysed 
in [l[ and also in 0] , where calls are not routed on direct links at all. 

Theorem 11.11 holds also in the case above where direct links are not used 
(i.e., each arriving call is allocated to the best among d indirect routes), with 
a suitably modified function F in the differential equation (|1.3p . Indeed, for 
< < C, we take instead 

Fk[i) = \9k-m - A<?fc(o - km + {k + i)i{k + 1), 

where the functions gk{£,) are amended by dropping the factor ^(C); -Fo(0 
and Fc{^) are modified in the same way. The proof is essentially identical, 
indeed slightly simpler in a few places. 

In (as well as in [8,] for a corresponding discrete-time model), the class 
of routing strategies choosing a path for a new call from among d random 
alternatives is called the GDAR (General Dynamic Alternative Routing) 
Algorithm. The particular model we study in this paper is called the BDAR 
(Balanced Dynamic Alternative Routing) Algorithm. The FDAR (First 
Dynamic Alternative Routing) Algorithm always chooses the first possible 
alternative two-link route among the d chosen. As in other models of this 
type, the 'power of two choices' phenomenon has been observed, that is, 
with the BDAR algorithm for d > 2, the capacity C = C{n) required to 
ensure that most calls are routed successfully is much smaller than with the 
FDAR algorithm. This phenomenon is not exhibited by the model studied 
in this paper, as proved in Theorem 1.1 of [sl], but it does occur in the variant 
discussed above where direct links are not used. 

In particular, see [t^, for the variant where there is capacity division into 
three parts, but with zero capacity for direct calls, the following is true after 
a 'burn-in' period of length O(logn). Suppose we use the FDAR algorithm 
and each indirect link has capacity C{n) ~ Q; io^giogn ' where a > is a 
constant. If q > 2/d, then there exists a constant K > Q such that the 
mean number of calls lost in an interval of length is o(l). If a < 2/d 
then for each K there exists a constant c > such that the mean number 
of calls lost in an interval of length is at least n'^. On the other hand, 
suppose we use the BDAR algorithm with d > 2 choices, and let K > be 
a constant. There exists a constant c > such that, if C(n) > '°^^g | " + c 
then the expected number of lost calls in an interval of length is o(l); 
and if C(n) < — c then the expected number of lost calls in such an 

interval is at least n^~^'^~°^^^ as n — t- oo. 

Our methods apply to any GDAR algorithm, and indeed any of the vari- 
ants discussed above. The law of large numbers for the BDAR algorithm 
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proved here is valid for the model without direct links in the parameter range 
considered in [7], i.e., with constant A and d, and C = C{n) = O(loglogn) 
for d>2, and C = C{n) = O (log n/ log log n) for d = 1. 

In [2], we carry out an improved analysis of the coupling introduced in 
this paper, to prove rapid convergence to equilibrium, and concentration of 
measure in equilibrium, for small enough arrival rate A. Also, for such A, we 
are able to show that ;^/^,a: is well-approximated by the unique fixed point 
of the differential equation (11. Sh . For d = 1, this proves the approximation 
suggested in [3], and it also provides an alternative proof of some of the 
results in 0], for small enough A. 

The rest of the paper is organised as follows. In Section 2 we develop a 
concentration of measure inequality that will be a fundamental ingredient of 
our proofs. In Section 3, we formally write down the generator of the Markov 
chain in question, and give an informal explanation of our proof strategy. In 
Section 4, we describe a simple coupling between two copies of the Markov 
chain, and show that, under the coupling, they do not get much further 
apart over time, according to suitable notions of distance. In Section 5, we 
use the coupling to show that nice functions of the Markov chain are well 
concentrated around their expectations. Section 6 contains estimates of the 
expectation of the generator of the Markov chain. In Section 7, we prove 
our main results. Theorems 11.11 and II. 2i In Section 8, we discuss the issue 
of initial conditions needed for Theorem 11.11 and 11.21 to apply. In Section 9 
we discuss how our techniques can be extended to analyse other models. 



2. Concentration inequalities 

We present some concentration of measure inequalities that will be used 
in our proofs, and may also be useful in the analysis of other Markov chains 
with similar properties. These inequalities generalise results presented in 0. 

Let X = {Xt)t^i+ be a discrete-time Markov chain with a discrete state 
space S and transition probabilities P{x, y) for x,y G S, where YlyeS v) — 
1 for each x £ S. We allow X to be lazy, that is we allow P{x, x) > for 
some X £ S. We assume that, for each x £ S, \{y : P{x,y) > 0}| < oo. 

This setting is natural, and many models in applied probability and com- 
binatorics fit into this framework, including those discussed in Section [H 

Let 

Q = = {oj = {ujQ,uji, . . .) : LOi £ S Vi}. 

Members a; of will correspond to possible paths of the chain X, in that 
Xi{u}) = uJi for i £ Then for each t £ Z+, Xt may be viewed as a 

random variable on a measurable space (fi, where J- = <T(U^Qj-i) and 
J-t = (T{Xi : i < t). The a- fields J-'t form the natural filtration for X. 

Let V{S) be the power set of the discrete set S. The law of the Markov 
chain can be interpreted as a probability measure P on and is de- 

termined uniquely by the transition matrix P together with the initial state 
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Xo = zo, according to 



P({w : ujj = Xj for all j < i}) = P(xj,Xj+i), 

i=o 

for each E S, for each z G Z"*". To be precise, this defines the 

law of (Xt) conditional on Xq = xq, and will be denoted by Fxq in what 
follows. Let P^{x,y) be the t-step transition probability from x to y, given 
inductively by 

P*(x,y) = ^P*-i(x,z)P(z,2/). 

Then, for ^ C 5, 



{Xt e A) 



ZO 



Let Ex^o denote the expectation operator corresponding to P^^f,; then E^.^, [/(Xj)] 
is the expectation of the function / with respect to measure d^^P^ ■ 



For t G Z"*" and / : 5" — )• M, define the function P*/ by 

(p'f){x) = Y,P\^.y)f{y). ^^s. 



Thus {P^f){x) = E^[/(Xt)] for each x. Also, for x G 5, let iV(2;) = {y g 5 : 
P(x,y) >0}. 

The following concentration of measure result for real-valued functions 
of Xt is presented to set the scene, and because it may prove to be of 
independent interest. We will not use it in the proof of Theorems ll.ll and ll.2[ 



thm. conca-general Theorem 2.1. Let P he the transition matrix of a discrete-time Markov 
chain with discrete state space S. Let f : S be a function. 

(i) Let {ai)i^z+ ^ sequence of positive constants such that for alii G Z, 



sup \Ex[fiXi)]-Ey[f{Xi)]\<a^. 

xeS,yeN{x) 
Then for all a > 0, xq £ S , and t > 0, 

F,,i\fiXt)-E,,[fiXt)]\ >a)< 2e-'^'/2(Eto^°?). 

(ii) More generally, let Sq be a non-empty subset of S, and let 
i^i)iez+ ^6 a sequence of positive constants such that, for all i G Z, 



(2.1) cond-gen 



(2.2) ineq. conca-general 



sup \Ex[f{X,i)]-Ey[f(Xi 

x,y&So.y&N(x) 



< Oii. 



(2.3) cond-gen-1 



Let Sq = {x £ So : y £ So whenever P{x,y) > 0}. Then for all xq G Sq, 
a > and t > 0, 

({|/(Xi)-E,„[/(Xi)]| > a}n{X, G 5° : 0<s< t-1}) < 2e-"'/2(Ea 

(2.4) 
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Our proof of Theorem l2.1l makes use of a concentration inequality from [13] • 
Let {Cl, T ^ P) be a probability space, with finite. (The arguments used here 
could be extended to many cases where Vl is countably infinite.) Let Q T 
be a (T-field of subsets of Then there exist disjoint sets Gi, . . . , Gm such 
that ri = U^j^Gr- and every set in Q can be written as a union of some of the 
sets Gr- Given a bounded random variable Z on (0, J^, P), the conditional 
supremum sup(Z \Q) oi Z m. Q is the ^-measurable function given by 

sup(Z I G){oj) = min _maxZ(£l'') = max Z{Cj'), (2.5) 



max Z{Cj), 

Cj'eGr 



where w € Gr- Thus sup(Z | Q) takes the value at u) equal to the maximum 
value of Z over the event Gr in Q containing a). 

The conditional range ran(Z) of Z in Q is the ^-measurable function 

ran(Z | Q) = sup(Z | G) + sup(-Z | G), (2.6) 

that is 



ran(Z | ^)(a)) 



max_ |2'(lji) — Z(a;2) 

(1)1, [1)2 SG,. 



where u> E Gr- Note that sup(Z|^) and ran{Z \ Q) are ^-measurable. 

Let t S N, let {0, fi} = J"o ^ -^i ^ • • • ^ -^t be a filtration in J^, and let 
Zq, - - - ,Zt he the martingale defined by Zi = E(Z| J^j) for each i = 0, . . . ,t- 
For each i, let rauj denote ran(Zj|J^j_i); by definition, ran, is an J-'i^i- 
measurable function. For each j, let the sum of squared conditional ranges 
be the random variable X]i=i ^an?, and let the maximum sum of squared 
conditional ranges be the supremum of the random variable that is 



f2 
' 3 



sup i?|(a}). 



The following result is Theorem 3.14 in 

Lemma 2.2. Lei Z be a bounded random variable on a finite probability 
space (0,-F,P) with E{Z) = m. Let {(Il,n} = <^ Ti C . . . C Tt be a 
filtration in T , and assume that Z is Tt-measurable. Then for any a > 0, 

P(|Z-m| > a) < 2e-2'^'/^?. 

More generally, for any a > and any value r|, 

P({|Z - m| > a} n < r1}) < 26-^'^"/'^?. 



Proof of Theorem [27TJ Let us start with (i). Let / : S — )• M be a function. 
Fix a time t G N, and an initial state xq G 5" and consider the evolution of 
Xt conditional on Xq = xq for t steps, that is until time t- Since we have 
assumed that there are only a finite number of possible transitions from any 
given X G S", we can build this process until time t on a finite probability 
space (fi, P^q): we can take O to be the finite set of all possible paths 
Xq, . . . .,xt of the process starting at time in state xq until time t- 
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thm. concb-general 



For each time j = 0, . . . ,t, let J^j = a{XQ, . . . ,Xj), the cr-field generated 
by Xq, . . . ,Xj; so J^o = {0, Cl}. Also, we let T = Pf 

Consider the random variable Z = f{Xt) : — t- M; note that f{Xt) is 
J"t-measurable. Also, for j = 0, . . . , t let Zj be given by 



Zi 



.J - E,,[f{Xt)\fj] = E,,[f{Xt)\Xo, ...,Xj] = iP'-^f){X,), 

where we have used the Markov property in the last equality. 

Fix 1 < J < t; we want to upper bound ran-,- = ran{Zj \ Tj-i). The 
(T- field J'j^i can be decomposed into events {Cj : oji = Xi for alH < i — 1}, 
for different possible paths xq, xi, . . . , Xj-i of X. Fix xi, . . . , Xj-i G S, and 
for X G 5 such that P(xj_i,x) > consider 



h{x) 



E,„[/(Xt)|Xo 



,X,=x] = {P'-^f){x). 



Note that Zj{Cj) G {h{x) : P{xj^i,x) > 0} for u such that Xj-i{u) 
It follows from (12.11) that, for such a). 



Xj-l- 



ranj{u) 



< 



sup \h{x) 

x,y:P{xj-i,x)>0,P{xj-i,y)>0 

2 sup |(P*-V)(x,_i)- 

x:P(xj — i,x)>0 



Hy)\ 



< 2a 



t-j- 



It follows that 



t-i 



i=0 



uniformly over oj €z Q. Part (i) of Theorem 12.11 now follows from Lemma [2.2[ 
To prove (ii), observe that the bound 

ranj(a}) = ran{Zj \ Fj-i){Cj) < 2at-j 

still holds on the event At = {u : Xj(oj) G Sq for j = 0, . . . , t — 1}. □ 

The next, more refined, concentration of measure result, is the one that 
we actually use in our proofs. 

Theorem 2.3. Let P be the transition matrix of a discrete-time Markov 
chain with discrete state space S. Let f : S be a function. Suppose the 
set Sq and numbers ai{x,y) (x,y €z Sq) satisfy the following conditions: for 
all i G and all x, y G Sq, 



\E,[f{X,)]-Ey[f{Xi)]\<a^{x,y). 

- {x £ Sq : y £ Sq whenever P{x,y) > 0}. For x G S and i G 
— 7- M &e given by 



Let S^ 

ax,i{y) = ai{x,y). 
Assume that, for some sequence {oii)i^z+ of positive constants, 

sup (Pal i){x) < al 
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thm.mart-b 



Let t > 0, and let (3 = 2 X]i=o ^i- Suppose also that a is such that 
2 sup sup ai{x,y)<a. 

0<i<t-l x£S^,y£N{x) 
Finally, let At = {u : Xg^uj) £ Sq : < s < t — 1} . Then, for all a > 0, 



(2.9) cond-gen-5 



{\fiXt)-E,,[fiXt)]\>a}nAt) <2e 



-aV(2/3{l+(aa/3;9)) 



(2.10) 



To prove Theorem 12.31 we use another result from 10(]. 

As before, let Z be a bounded random variable on a finite probability 
space (i^,J",P). Fix t G N, and let {0, (7} = J"o ^ • • • ^ C J^- be a 
filtration. For i = 0, ... ,t, let Zi = K{Z \ Ti). For i = 1, . . . , t, define the 
conditional variance 

raviZi I Ti-i) = tdZi - E{Zi I J-i_i))2 I 



varj 

Further, let V = Yll=i '^^^i: the sum of conditional variances. Also, for each 
such i = 1, . . . , t, define the i-th conditional deviation 

devj = sup(|Zi - I Ti^i), 

and let the conditional deviation be dev = maxjdevj. Note that V and 
dev are random variables in (Q,J^,¥). The following result is a 'two-sided' 
version (and a simple consequence) of Theorem 3.15 in [lo| . 

Lemma 2.4. Let Z be a random variable on a finite probability space (il, J-, P) 
with E(Z) = m. Let {0, il} = /"o ^ -^i ^ • • • ^ -^t be a filtration in T , and 
assume that Z is J^t -measurable. Let a = max^^^ dev(a;), the maximum 



conditional deviation. Let (3 



max- 



V{Cj), the maximum sum of condi- 



tional variances. Assume that a and (3 are finite. Then for any a> 0, 
F{\Z -m\>a)< 2e-"V(2/3(i+(W3/3)), 
More generally, for any a > and any values a, (3 > 0, 

P({|Z -m\>a}n{V</3}n {dev < a}) < 2e-"'/(2/3(i+(W3/3)) _ 

Proof of Theorem 12.31 We start, as in the proof of the previous theorem, 
by assuming that 5*0 = S. Let / : S" — t- R be a function. Fix a time 
t G N, and an xq G 5; consider the evolution of X = {Xt)t>o conditional on 
Xq = xq for t steps, that is until time t. Again this process can be supported 
by a finite probability space 

For each j = 0, . . . ,t let j^j = cr{XQ, . . . , Xj), the cr-field generated by 
Xq, . . . ,Xj; so To = {0,17}. Also, let T = Tt- We consider the random 
variable Z = f{Xt) : O — )• M. And, for j = 0, . . . , t, Zj is given by 

Z, = KAf{Xt)\Pj\ = ]E,o[/(Xi)|Xo, ...,x,] = {P'-^f){X,). 

We want to apply Lemma 12.41 so we need to calculate the conditional vari- 
ances varj. We use the fact that the variance of a random variable Y is 
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equal to ^K{Y — Y)"^ , where Y is another random variable with the same 
distribution as Y and independent of Y. 

Fix I < j < t and xi, . . . , xj-i £ S, and for x G consider 

h{x)=E^,[f{Xt)\Xo=xo,...,Xj^i = xj^i,Xj =x] = iP''^f){x). 

Then, for a) such that = Zj{uj) £ {h{x) : P{xj-i,x) > 0}, 

so that 

vaijiLj) = ^^P{xj-i,x)P{xj-i,y){h{x) - h{y)f 

x,y 

x,y 
x,y 

+ ^P(x,_i,x)P(x,_i,y)((P*-V)(^i-i) - 
x,y 

< 2 Yl P{x,-i,x){{P'-^f){x)-{P'-^f){x,^r)f 

x:P(xj^i,x)>0 

X 

by assumption ()2.8p . It follows that 

t-i 

/3</3 = 2^a2. 

We now bound dev = maxi<j<j devj. For u such that = Xj-i, 

dev.iCj) = sup 

x:P{xj-i,x)>0 

sup |(P*-^/)(x) - 

x:P{xj-i,x)>0 

< sup - {P'-'f){y)\ 

x:P{xj -l,x)>0 y 

< sup j;P(x,_i,y)|(P*-V)(x) - (P*--'/)(x,_i)| 

a;:P(xj_i,a;)>0 

+ Y,P{x,.,,y)\iP'-^f){y) - (P*-V)(x,-i)| 
y 

< 2 sup |(P*-V)(a;)-(P*-^/)(x,-i)| <a, 

a;:P{a;'j_i,a;)>0 

for each 1 < j < t, by assumption ()2.9p . Therefore a < a. 
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Theorem 12.31 now follows from the first statement in Lemma 12.41 in the 
case where 5*0 = S. In general, the above bounds on V and dev hold on the 
event At = {uj : Xi{uj) £ Sq for i = 0, . . . ,t — 1} , and so the full statement 
of Theorem 12.31 also follows from the second inequality in Lemma 12.41 □ 

We now prove that the expectation of a well concentrated function / 
multiplied by an indicator function approximately factorises, with bounds 
in terms of bounds on / and its deviations from its mean. This result will 
be used several times in our proof of Theorems 11.11 and 11.21 

lem.expec Lemma 2.5. Let {Q,J^,¥) be a probability space and let X be a random 
variable on (il, P) taking values in a measurable space {S,S). Let f : S ^ 
M. be a measurable function, and suppose that ¥{\f{X) — E f{X)\ > a) < b 
and P(|/(X)| < c) = 1. Let AgS. Then 

E[IxeAf{X)]-F{X G A)E[f{X)]\ < aP(X e A)+bc. 

Proof. We have 

E[IxeAf{X)] = E[IxeA/(^)II|/(x)-E[/(x)]|<a] 

+ E[IxeAf{X)I\f(^x)-E[f{X)]\>a\- 

Now, 

E[IxeAf{X)I\f^x)-nf{x)]\<a] < E[IxeA{nf{X)]+a)] 

= F{X G A) E[f{X)] + F{X G A)a, 

and 

E[IxeAf{X)I\f^xyEmx)]\<a] > E[IxeA(E[/(X)] - a)] 

= F{X e A) E[f{X)] - F{X G A)a. 

Also, E[IxeA/(^)I|/(x)-E[/{x)]|>a] < cE[IxeAl|/(x)-E[/(x)]|>a] < cb and 

^[^XeAfiX)I\f(^X)-E[f{X)]\>a] > -cE[IxeA^\f(X)-E[f{X)]\>a] > -cb- 

The result follows. □ 



S: gen 



3. Generator of the Markov process 

We now return to the routing model described in the introduction. 

Recall that fv,j{x) denotes the number of links with one end v carrying 
exactly j calls, and that luv{x) = 1 if x{{u, v}) = j and 1^(3;) = otherwise, 
for all u,v,j. Also, we let I^i{x) = 1 if x{{u,v}) < j and I^i{x) = 
otherwise. Further, we define = 1 if x{{u,w}) V x{{v,w}) < j, and 

^uv,w{x) = otherwise. 

Let A be the generator operator of the Markov process X. By standard 
theory of Markov chains, for each t > 0, each v £ Vn and each j G {0, . . . , C}, 

dE[f,j{Xt)] 



dt 



:E[Af,jiXt) 
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SO in order to prove Theorem 11.11 we will need to be able to approximate 

For X £ S and < j < C, we can write 

Afv,j{x) = \fvj-i{x) - Xfv,j{x) + Xgvj~i{x) - Xgv,j{x) 

- jfv,j{x) + {j + l)fvj + l{x), 

Afv,o{x) = -Xfv,o{^) - Xg^fi{x) + fv,i{x), 
Afv,c{x) = Xfy^c-iix) + Xg-a^c-iix) - Cfv^c{x), 

where the g^j (x) represent contributions due to alternatively routed arrivals 
with one end v, and are given, for j = 0, . . . , C — 1, by 

^ ' r = l MjW S = \ s=T- + l 

d C-\ r-l d 

r=l lijW j=j-|-l s=l s=r+l 

r=l n,i;',Wr s=\ s=r+l 

+ E E ^^nA E ^f^'/n(i - n(i-i.-:i)]Ai) 

r=l u,ii',Wr i=j+l s=l s=r+l 

with W = (Wi, . . . j-Wd) G = (Wl, . . . ,10^-1, Wr+1, • • • ilL'rf) G K''^^ 

Here, ^ denotes the sum over all n / and over all . . . ^w^i such 
that each Wr / n, t>, and ^„ ^/ denotes the sum over all u ^ v, v' ^ u,v 
and over all wi, . . . , Wr-i-,Wr+i-, . . . ,Wd such that each wj ^ u, v' . 

In ()3.ip . the first term is the probability that the direct link chosen for a 
new call with one end v is blocked and, on the two-link path selected for the 
call, the link including v has load j, while its partner link has load at most 
j. The second term is the probability that the direct link chosen for a new 
call with one end v is blocked and, on the two-link path selected for the call, 
the link including v has load j and its partner link has load greater than j. 
The third term is the probability that v is chosen as an intermediate node 
for a call blocked from its direct link, the route through v is the best out of 
the d routes selected and j is the maximum load of a link on the route. The 
fourth term is the probability that v is chosen as an intermediate node for 
a call blocked from its direct link, the route through v is the best of the d 
routes selected, and j is not the maximum load of a link on the route. 

In particular, when d = 1, for j = 0, . . . , C — 1, for each v G Vn, 

— ^ \^ (]\'^ P T[<C-l^TfC'p IT^C'-l^ 

Tl ^ 

U,WGVn 
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Furthermore, when d = 2, then g^j is a sum of 8 contributions. These 
contributions correspond to the case where v is an end node and the case 
where v is an intermediate node for a call. In the case where v is an end 
node, we have a subcase where the first route of those selected is allocated 
to a new call and a subcase where the second route of those selected is 
allocated to a new call. We further need to distinguish a subcase where 
the link ending in v has the maximum load, and a subcase where the link 
ending in v does not have the maximum load on the route allocated to a new 
call. In the case where v is an intermediate node, we need to distinguish a 
subcase where v is the first alternative node selected and a subcase where v 
is the second alternative node selected. Also, we have a subcase where the 
link with load j has the maximum load, and a subcase where a link with 
load j does not have the maximum load on the route allocated to a new call. 

The contribution due to the case where v is an end node, the first route 
of those selected is allocated to a new call, and the link ending in v has the 
maximum load is of the form 

— - — V 

U,W\,W2&V„ 

The contribution due to the case where v is an end node, the second route 
of those selected is allocated to a new call, and the link ending in v does not 
have the maximum load is of the form 

C-l 

\ ^ ttC ttj \ ^ iri (-1 •[r<i \ 

in — 2)"^ ^ ■"■«1'^1'UI2 / J '^UW2\ '^UV,Wl)- 

U,W-i_,W2eVn i=3 + '^ 

The contribution due to the case where v is an intermediate node and is 
selected first, and where a link with load j has the maximum load on the 
route allocated to a new call is of the form 

The remaining contributions can be expressed analogously, and the form of 
g for (i > 2 is derived similarly. 

We note that each g^^j is a sum of products of indicators of sets of load 
vectors with properties pertaining to loads of particular links. Our plan is 
to justify the intiuition that, subject to suitable initial conditions, the loads 
on different links behave nearly as iid random variables at each time i, and 
the precise estimates we use involve sums over reasonably large collections 
of links. We need several specific manifestations of this near-independence 
and symmetry. First, the geometry of the network is not important; this 
means that, for fixed nodes u and the loads on links uw and vw are not 
strongly correlated, on average over u;, so that the average value of '^uw^w 
over w is close to the product of the average values of I:^^ and I^^. (In 
other words, the function (^"^ defined earlier is small.) Secondly, for fixed 
nodes u and f , the loads on links incident on u have approximately the 
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same distribution as the loads on links incident on v. (This means that 
the function 0^ small.) Thirdly, we require that the alternatively routed 
calls are fairly uniformly distributed over the network. (This implies that 
the function (j)^ is small.) Finally, we will show that each fyj{Xt) is well 
concentrated around its expectation, which then implies that, for fixed nodes 
u and V, E[fu,j{Xt)fv,j{Xt)] is approximately equal to E[fujiXt)] E[fvj{Xt)]. 
Naturally, quantitative versions of these properties will need to be assumed 
to hold at time 0, and we will show that they are maintained throughout 
the time period of interest. This will then allow us to express E[(7^j(Xt)] 
as a (scaled) sum of products of terms of the form :f^fv,i{t), and hence 
lead to approximate differential equations satisfied by the K[fvj{Xt)] for 
J = 0, 1, . . . , C, for each v £ Vn, expressed in terms of themselves. 

Let fv,<j{x) = J2i<j fv,ii^)- Let fvj{t) = E[/„j(Xj)], and let fv,<j{t) = 
E[h,<j{Xt)]. Let lL{t) = E[lL(Xt)] and let l|j(t) = ¥.[lti{Xt)]. We will 
show that the expectation of the first term in (13. ip with respect to the law 
of Xt is approximately 

^ ' r=l u w s=l 

d 

X n {^-^tt^\)){xt) 

s=r+l 

r-1 



^ E E ^ [( E 11(1 - i^^^jf^: 



^ ' r=l u 



d 

X 

s=r+l 



=r+l 

^ Tr^df'',c{t)u,{t)f,,<,{t) E(i - {fv,<m'y-' 

^ ' r=l 

x(l-(A<.-iW)Y-^ 
Handling the remaining terms in a similar way, we will prove that 



2 



n9vAXt)] « , _ ^s2d Ag(t)Aj(t)/.,<j(t) E(i - {fv,<3{t)?y-^ 



x(i-(^<,„i(t))Y- 

2 



c-i d 

+ (^_i)2d f^,c{t)f.At) E fvAt)Y.(^ - ifvMt)fY-' 



x(i-(A,<,_i(t))^) 
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ffv,cit)fv,j{t)fy^<j{t) 



sec : couple 



^, (i-(/.,<.-i(t))Y-(i-(/.,<,(0)Y 

f.,<j{t)r - (.f.,<j-i(.t)r 

2 

^ ' i=j+l 

(i-(/.,<.^i(t))y-(i-(/.,<.(t))y 
f.,<iit)r - ifv,<i-i{tw 

Hence we will see that the functions K[f,uj{Xt)] (j = 0, . . . , C, w G Vn) 
approximately solve the differential equation (II. 3p . As the fvj{Xt) are well 
concentrated around their expectations, Theorems II . II and II . 2 1 will follow. 

4. Coupling 

In this section, we describe and analyse a natural coupling between two 
copies of process X. We start by defining notions of 'distance' between the 
two copies, with the aim of showing that the expected distance grows slowly 
in time, at least for time 0(1). 

Given two load vectors x, y, the ^i-distance between them is 

11^ ~ = E ~ + E ^(e,'^) - y(e,?i')|, 

which measures the sum of the differences between x and y in loads of all 
the different routes. Then || • || is a metric on S. 
For V £ Vn we will also consider functions 

= '^\x{e,0) -y{e,0)\+ '^'^\x{e,w) -y{e,w)\ 

e:v£e ewGew^e 



+ E \x{{u,w},v) -y{{u,w},v)\. 

Then 2||rE — gives an upper bound on the sum of the differences between 
the loads of links {v,w}, w v (i.e. links around node v) in x and y. 

For the remainder of this section, and subsequently, we will work with a 
jump chain X corresponding to X. The chain we use is not the standard 
embedded chain but a slower moving version that will often not change its 
state at a given step. Since we only work on the set S, we can use a jump 
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chain with transitions defined as follows. Given that the current state, at 
time t G Z+, is x £ S", the next event is an arrival with probability 

= A© + L6A(;)J 

and a potential departure with probability l—p{X, n). Given that the event is 
an arrival, each pair of endpoints u,v is chosen with probability 1/(2)) then 
each d-tuple of intermediate nodes is chosen with probability (n — 2)~'^, and 
the call is routed along the two-link route chosen first among the d selected 
that minimises the maximum load of a link. Given that the event is a 
potential departure, the calls currently in the system are enumerated from 
1 up to at most [6A (2) J , and then a number is chosen uniformly at random 
from {1, . . . , [6A(2)J}. If there is a call assigned to this number, it departs, 
and otherwise nothing happens. 

Let 5*0 C 5 be the set of states x such that ||x||i <^\{^. Recall also that 
Si C 5*0 is the set of states x such that ||x||i < 2A(2). We will be interested 
in the evolution of the chain starting from S\ and before it leaves Sq. 

Consider the following family of Markovian couplings (X^" , Y^^ ) of pairs 
of copies X^'",y^'' of the discrete jump chain starting from states xo,yo re- 
spectively, where xq, yo £ Sq. (In what follows, we will drop the superscripts 
xo, yo from the notation and refer simply to X and Y .) 

Let t > 0, and let x. y be both in S. Given that Xt-i = x and Yj-i = y, 
the transition at time t (from state {Xt-i,Yt-i) to {Xt,Yt)) is an arrival 
in both X and 1", or a potential departure in both X and Y . Given that 
the transition is an arrival, we choose the same call endpoints and the same 
d-tuple of intermediate nodes in both. Also, given that the transition is a po- 
tential departure, we pair the calls occupying the same route in both X and 
F, as much as possible. We also pair off the remaining calls arbitrarily, as 
much as possible. (We can pair off all the calls if||a;||i = ||y||i, but otherwise 
some calls will remain unpaired in the process that currently has more calls.) 
Then we always choose the paired calls for simultaneous departure; any calls 
that are unpaired will depart alone. This can be achieved by assigning to 
each pair, and to each unpaired call, a distinct number in {1, . . . , [6A(2)J}- 
If the transition at time t is a potential departure, we choose the same uni- 
formly random number from {!,..., [6A(2)J} for both X and Y. If the 
number corresponds to a pair of calls, both depart; if it corresponds to an 
unpaired call, this call departs; otherwise, nothing happens. 

Suppose now Xt~i = ^ and Yf-i = y where at least one of x and y, say 
X, is not in S. We then let X^ = X^^i = x, while Yt is obtained from Yf-i 
by following the transition probabilities of the jump chain from y if y G 5 
OTYt = Yt.i = yiiy^S. 

The process W = (Wt) given by Wt = {Xt, ^t) is a Markov chain adapted 
to its natural filtration, Qt = cr{Xs, Yg : s < t). 
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Note that, on the event that the jump at time t is a potential departure, 
^ ll^~y||i- (The distance remains unchanged if paired calls from 
the same route depart or if there is no departure at all; it decreases by 2 if 
paired calls on different routes depart, and decreases by 1 if an unpaired call 
departs.) Furthermore, the distance between X and Y can only increase by 
2 at an event, and this can only happen if the jump is an arrival and if we 
select at least one of the links where Xt-i and Yt-i differ. This happens 
with probability at most 



2d+l 
1^ 



3d 



\Xt^i{e)-Yt-,{e) < THY E \Xt^i{e)-Yt,i{e) 



and Y.eeE„ \^t-i{e) - Yt-i{e)\ is equal to 

^ \{Xt-i{{u,v],^)-Yt-i{{u,v],{))) 

+ Y{Xt-i{{v,w],u)-Yt-i{{v,w],u)) 

+ i^t-i{{u,w},v) -Yt-i{{u,w},v)) 



< 2\\Xt^ 



Yt-i\\i. 



It follows that, uniformly over all x,y £ S, 



E 



\Xt-Yt\\i I Xt-i = x,Yt-i = y 



< 1 + 



I2d\ 



X - y\\i- 



We have assumed that Xq = xq and Yq = yo, where xo,yo £ Sq, that is 
ll^olli < 4A(2) and ||yo||i < "i^^Q- Note that, whenever > 4.X{'^), 

provided n > max (3, j), 

AQ) . 1 



t\\i 



Xt- 



and 



t\\i 



IX- 



t-1 111 



nil 



1 I X_i] < 



-1 I X-i] > 



< 



AQ) + L6AQ)J - 6 
4Am 



4 1 
> - > -. 



A(2) + L6A(^)J -7-2- 

Therefore, by standard inequalities (see, for instance. Lemma 2.5 in [7], with 
p = 1/6, q = 1/2 and a = 2A(2)), for any constant c > 0, 



3t < cn^ : \\Xt\\i V ||yf||i > 6A 



n 



< 2cn^ - 
3 



Let Ag be the event {X^j £ S,Yu ^ S for all u < s}. Then, for xq, yo £ Sq, 

E(||X-yt||ilA,_i) = E[K{\\Xt-Yt\\ilA,.,\Qt-i) 



1 + 
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< (i + |j^)e(||X,_i-F*_i||iI^,_,). 



By induction, 



K{\\Xt - YtWiU,.,) < (l + ^) Vo - 2/o||i 



Since the chain stops once it leaves S and each jump changes ||^t||i and 
||yt||i by at most 1, on the event At-i, \\Xt - Yt\\i < 2(6A(2) + 1). Hence, 
for any constant c < n, and any t < cn?, 

E{\\Xt-Yt\\i) = E{\\Xt-Yt\\ilA,_,) + E{\\Xt-Yt\\ii^-^A,.,)) 
< (l + |f)*IN-..ll. + 2(6AQ+l)2c„= 

(^ + "(fyj lko-2/o||i + 14An5(^-J 

If n > max(1000, 1/A), this last term is at most 1. Therefore, for n > 
max(1000, 1/A, c), and t < an?, 



E(||Xi-yt||i)<2(l + ^V 



(2) 



\xq - yolli, 



(4.1) I eq.dist 



uniformly over all possible starting states xo,yo in "So- 
Let u be a vertex. While in S, we can only change the loads of links at v 
(i.e. links {v, w}, for w ^ v) \i we choose a link with end v at a jump time. 
Also, we can only make — bigger than — at an arrival 

time, if either we pick one of the links {f , (if any) in which Xt^i and It-i 
differ, or if we pick a link {u, w} (where u ^ v and w ^ u,v) in which 
and Yt_i differ, and also vertex v as an endpoint or an intermediate node 
for a new call. The former happens with conditional probability at most 

'^l±lY.\^t-i{{vM)-yt-i{{vM) 

V2/ w^v 



\2/ w^v 



{{v,w}) 



and 



J2 I {Xt-i{{v, w},0)- Yt^ii{v, w}, 0)) 



X] (.^t-i{{u,v},w) -Yt-i{{u,v},w)) 
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+ J2 (^t-^(i'^M,v)-Yt^i{{u,w},v) <2||lt_i-yt_i||„. 

The latter happens with conditional probability at most 
1 2(f 



J2 \xt^i{{u,w})-yt-i{{u,w})\ 



provided n > 4. Also, always with probability 1, 

\\Xt-Yt\\v<\\Xt-i-Yt_i\\v + 2. 

Then, for n > 4, 

Ud 



E{\\Xt-Yt\UlA,_,) < {l + -^)E{\\Xt-i-Yt-i\UlA,_, 



and so 

K{\\Xt-YtUA,.r) < (l + ^)Vo-yo||. 

+ ^?jr(l + 7Hr) lko-yo||i- 

Hence, as before, for all n > max(1000, 1/A, c) and t < cn'^, for each 
V G Vn, and for all xq, yo € ^o. 



+ -pr(l + 7m) lko-2/olli. (4.2) 



"'V2; ' V2, 

Recall that, for a load vector x, fv,k{x) is the number of links around v 
carrying exactly k calls. Similarly, fv,<k{x) = Ylii<k fv,i{x) is the number of 
links {v, w}, w ^ V, such that x{{v, w}) < k; that is, the number of links 
around v carrying at most k calls. Let P denote the transition matrix of 
the jump chain (Xt) restricted to S. 

Note that, for each v, k, and each x,y e S, 

\fv,k{^) - fvAy)\ ^ ^\x{{v,w})-y{{v,w})\<2\\x-y\\y, 

\fv,<k{x) - fv,<kiy)\ < ^\x{{v,w}) -y{{v,w})\ <2\\x-y\\y. 
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Hence, for xq, yo G ^o, for n > max(1000, 1/A, c), for t < cn^ and each v, k, 
\iP'fv,k){xo) - iP'fv,k){yo)\ < E \f,4Xt) f,,k{Yt)\ < 2E \\Xt - YtWv 

^U^t ( 12d\*-i, „ 
+ 1+ ll^o-yolli. (4.3) 



Similarly, for xo,yo G Sq, n > max(1000, 1/A, c), t < cr? and each 
|(P%<fe)M-(PV.,<ik)(yo)| < 4(l + i^)Vo-2/o||. 

No-»lli.(4.4) 



Given u^v and j, fc G {0, 1, . . . , C}, set 
(n-2) 



^ a r— i 

fu,v,.j,k = — nx^^_i ^Mw^ Kwi. tt xi^-'-^-'^^'i^s-'^'^^ 



s=r+l WstWg 

d r-1 



J a r— i 

fu,v,<j,k = / _ 2)2d-i yz X] ^"'i''- X] ^^'"r n ~ 



X 

s=r-+l Ws,w' 



where the sums are over all Wr,w'j.,Ws,w'g ^ u,v. Then 
1 

in -'2) 

d 



fu,v,j,k — _n\2d-l^^'^'i 'Hiv)Uv,k \v) 



X ^ ((n - 2)2 - - lU){fv,<j - l^i)f ' 

r=l 

X ((n - 2)2 - - i^r')(/.,<.-i - i^r'))'"', 

fu,v,<j,k = (jj, _ 2)2'^"^ ifu,<j ~ ^uv)ifv,k ~ ^uv) 

^ ((n - 2)2 - - 4^)(/.,<, - l^^)y~' 

r=l 

((n - 2)2 - - 4^-i)(/.,<,_i - l^r'))''^- 



X 
r-=l 



Malwina J. LUCZAK 



sec : conc-route 



In the case d = 1, we have 

fu,v,<j,k{x^ 



n-2 



so that 



fu,v,j,k 
fu,v,<j,k 



n-2 



ifu,j-iL)ifv,k-iL 



—:^{fu,<j - l^i){fv,k -\v)- 



n-2 



Then 



\fu,v,j,k{,x) fu,v,j,k(,y)\ 

^ - lUx))\fuA^) - iLix) - ifuAy) - lL{y))\ 



< 



n 



+ 



1 



{fuAv) - KM)\fv,k{^) - tv{x) - iUiy) - tv{y))\ 



n-2 

< 2\\x -y\\u + 2\\x -y\\v 
Similarly, 

\fu,v,<j,k{^) - fu,v,<j,k{^)\ ^ - y\\u + 2\\x - y\\y. 

A calculation similar to the one above shows that, for any d > 1, if / is 
one of the functions fu,v,j,k, fu,v,<j,k^ then 

\f{x)-f{y)\<2d\\\x-y\\u + \\x-y\U). 

As before, if / is one of the functions fu,v,j,k^ fu,v,<j,k, then, for all xq, yo £ 
Sq, all n > max(1000, 1/A, c), and all t < cn"^, 

\{P'f){xo)-{P'f){yo)\ 

w 



I2d\t 



< U^(l + ^\\\xQ-yQ\\y + Ad^{l + -r^ 



Fo - yolk 



+ 



l2SdH 



n 



1 + 



12(i\*-i 



(2) 



Fo - yolli- 



(4.5) 



5. Concentration of measure for the routing model 



We will now apply Theorem 12.31 to the jump Markov chain X and func- 
tions fy^k: fv,<ki fu,v,j,ky fu,v,<j,k' 

From now on, we assume that our process starts in some fixed state 
-^0 = Xq G 5i. We write P and E when discussing probabilities relating 
to X, instead of and ^xq, which was convenient in the derivation of the 
concentration inequalities in Section [2j 
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We start with the functions f^^k- By ()4.3p . for all x,y £ Sq, we can take 

ai(x,y) = 4(1 + — ^1 \\x-y\U + ^^[l + -T^] \\x - y\\i, 

for i < cn? and n > max(1000, 1/A,c). 

Note that, for all x,y such that P{x,y) > 0, we have \\x — y\\v < 1 and 
\\x — 2/11 1 < 1- Also, for each x G 5o, 

, 2 + d 

y;||a;-y||„>0 

It follows that for each x G So, for i < cr? and n > max(1000, 1/A, c), 

~ n \ n J ~ n 

So we can take = d^c+if ^gedc ^ ^^^^ so /3 < 2^8^4(0 + l)3ne96'^^ for 
t < cn^. Also, for n > max(1000, 1/A, 64d^c) and t < cn^, we can take 

n 

By Theorem 12.31 for ^ ^ cn^, 

{\fv,k{Xt) - nfv,k{Xt)]\ > a} n {X, G : < s < t - l}j 

< 2e~''^/(^^^'^''('^+^)^"''^'*''''+^^''^**''''"/^) (5 1) 

Similarly, 

P ({|/„,<fc(Xt) - E[/,,<fc(lt)]| > a} n {X, G 5g : < s < t - 1}) 

< 2e~''^/(^"'^^(''+^)^"'''"'''''+^^''''*''''''/^) (5 2) 
We now consider functions fu,v,j,k and fu,v,<j,k- 

By for all 2;,y G So, 

aiix,y) = M'^ (l + \\x - y\\^ + 4d^ (1 + \\x - y\\u 



\2) 

128(i^w I2d\i-^. 



.2; - (2)^ " 

for i < en? and n > max(1000, 1/A, c). This leads to 

~ n \ n J ~ n ' 

for each x G 5o, for t < cin? and n > max(1000, 1/A, 64(i^c). 
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So we can take af = 2"^'^'(c+i)% 96dc^ ^^^^ ^ < 2^^d^{c+lfne^^'^^. Also, 
for n > max(1000, 1/A, GAcPc), we can take 



n 



By Theorem 12.31 

F - > a} n {X, G 5° : < s < t - 1}^ 

< 2e"''^/(^^^'^*(''+^)^"^'"^'^''+^^'^^''''*''''"/^) (5 3) 

and, similarly, 

F {{\fu,v,<j,k{Xt) - E[fu,,,<j,km]\ > a} n {X, G S^ < s < t - 1}) 

< 2e~"^/(^^^'^*('^+^)^"^^'''^''+^'''^^'^''*'^''"/^) (5 4) 
>2AQ), 



Now, whenever 

F{\\Xt\\i - \\Xt-i\\i = l\Xt-i) < 

and 

P(||Xt||i-||Xj_i||i = -l|Xt_i)> 



< 



1 

6' 



A(^) + L6A(^)J 

2A(^) ^ 2 
A(^) + L6AQ)J -7- 



By, for instance. Lemma 2.5 in Q, it follows that, for each xq G 5*1, 

t < cnMt 5g) < cn2(7/12)2M2)-i < cn^e""/^ < e-"/^ (5.5) 



provided n > max(1000, 1/A, 64d^c). 

Hence, for n > max(1000, l/X,6Ad?c), t < cn?, and a < n, we have 



-n/8 



\f.,km-nfv,k{^t)]\>aj 

< 4g-aV220d4(c+l)3ne96d=_ 

Similarly, we have, for n > max(1000, 1/A, 64(i^c), t < cv? ^ and a < n, 

r{\f.MXt)-nf.MXt)]\>a) < 4e-V2"^^(^+^)^-""(5.7) 
P {\fu,v,j,k{^t) - E[/„,,,,- > «) < 4e-'^'/2'''^'(^+i)'-''"(5.8) 

We also have 

P (3?; G K, t < cn^. A: G {0, . . . , C7}, \fy,k{Xt) - E[f.,,k{Xt)]\ > a) 

< 4cCn3e-'^'/2'°'^'(^+l)'-^'''^ (5.10) 
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To relate together the continuous-time process X and the discrete-time 
chain X, note that departures in X can be represented by a Poisson process 
of potential departures with rate [6A(2)J, together with a suitable process 
of 'choices' defined as follows. Since X stops as soon as it exits while X 
is moving, at each potential departure we can number all the calls in the 
system from 1 up to at most [6A (2) J • We then draw a number uniformly 
from {!,..., [6A(2)J}; if the number selected corresponds to a call, that call 
departs, and otherwise nothing happens. For this representation the mean 
number of events in X during the interval [0, io] is 



2;' 

and the events correspond precisely to the jumps of X. Let = c/8A, so 
that the mean number of events in [0, to] is at most ^cn^. Then, by standard 
concentration inequalities for Poisson random variables, the probability that 
there are more than en? events of X during [0,to] is at most e"'^" 

By the above, if n > max(1000, 1/A, 512(i^ Ato, l/8Aio), then for a < n, 
and any t < to; 

P [\U,,{X,) - nfv,k{Xt)]\ >a) < 4e-V2-<^^(c-M)3ne-- + ^-cnV6 

< 4g-aV220(i4(8Ato + l)3neSOO''^*0 _^ ^-n/& 

< 8e-"'/2'°'^''(^^*o+^^''''=*°°''^*°. (5.11) 
Also, under the same conditions, we have 

P(l/.,<fe(X*)-E[A,<,(XO]|>a) <8e-«V2-ci^(8Ato+i)Wooc<-o^ (5_^2) |eq.conc-f7a 



- {\fu,.,,k{Xt)-nk,v,,,k{Xt)]\ >a)< 8e-«V2-cl«(8A*o+i)3neBoo..*o^ (5_^3^ |eq.conc-f8 
F {\fu,v,<j,k{Xt) - nfu,v,<j,k{Xt)]\ >a)< 8e-«V2-<i«(8Ato+i)^ 



•)3jjg800dAto 



(5.14) eq.conc-fSa 

The above inequalities are all uniform in cZ > 1, u,!; G V^, and j, G 
{0, 1, . . . , C}. In particular, if we set 

1 

^ ~ 224^8 (SAto + l)3e800dAto ' 

and take a = ^y^logn, then we obtain that, for n > max(1000, 1/A, hl2(fi\to, l/8Ato), 
t <tQ, d>\, and each n, v and j, k, 

p(|/,,fc(Xt)-E[/,,fc(X4)]| >iv^logn) < 8e-T'°s'".(5.15) 

F {\fv,<km - nfv,<kiXt)]\ > ^V^logn) < 8e-^^°e^".(5.i6) 

F {\fu,v,j,km - E[U,,j,k{Xt)]\ > ^V^logn) < 8e-^'°s'" (5.17) 
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lem.gen-expec 



{\fu,v,<jAXt)-nfu,v,<jAXt)]\>lV^\ogn^ < 8e-^l°g'".(5.18) 



Hence 

P (3v €Vn,t< cn\ A; G {0, . . . , C}, \h,k{Xt) - E[/,,fe(Xt)]| > ^v^logn) 

< 8cCn3e-T^°e'". (5.19) 

6. Expectation of the generator 

As before, we assume that our process starts in some fixed state xq & Si, 
and we consider the law of the process started in this state, and running 
until some time to > 0. 

Recall that 



224^8 (8Ato + l)3e800'^^*o 



and 



no(A, d, C, to) = max (i, W^X + 1)U\C + if 4, e^M^'^'M 

Note that no(A, d, C, to) ^ e^'^' ^ 1000, for any positive integers d and C 
and positive reals A and to, so that the bounds of the previous section hold 
for ah n > no(A, d, C, to)- We also have e-T'°s " < n ^, an inequality we 
shall use freely from now on. 

Recall the definitions of (p^, cj)'^ and (p^ from ()1.4p - ()1.6p . and that 4> = 
max{(/)^, (/)2, (/)3}. Set (j) = max{(/)^, </>2}. Recall also that 

d 

gjiO = 2acwn<j)Y.(^-a<j?Y~\i-a<j--^f 



\2\d-r 



r=l 
C-1 



\2\d-r 



+2mm E mY.(^-ci<irY-'{l-a<^-m" 

i=j+l r=l 

Our first aim in this section is to show that, provided E,(p{Xt) is small, 
K[gTjj{Xt)] is close to (n — l)gj{z^), where is the vector with components 
zt{v,j) = (n - 1)"^ E[fvj{Xt)], for j G {0, . . . , C}. We then go on to show 
that, if (p{xo) is small, then also K(j){Xt) is small for all t < to- 

Lemma 6.1. For all t < to, for each v €z Vn and each j £ {0, . . . , C}, 

\ng,,,{Xt)]-{n-l)g,{z^)\ 
< 22d\C + lfn(<P{xo) + ^1^^.208(A+i)d2(c+i)3to 



n 



provided n > no(A, d, C, to)> where z^ is the vector with components zt{v,j) 
(j G {0, . . . , C}), with zt{v,j) = (n - l)-i nfvAXt)]. 

If d = 1 and n > no(A, 1, C, to), we have the improved bound 

\ng,,,{Xt)]-in-l)g,{z^)\ 
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3„r^r^„^ , 31og^\ 208(A+l)(C+l)to 



< 22{C + iyn((l){xo) + ^^)e 



The lemma above will follow immediately from two other lemmas, the 
first of which is as follows. 



lem.gen-exp Lemma 6.2. For any n > no{X,d,C,tQ), v G Vn and j G {0, . . . ,C}, 

E[g„j{Xt)]-in-l)gj{z^) < 8d^{C + lfnE[4>{Xt)] + 16d\C + l)V^logn. 

Proof. Suppose that n > no(A, d, C, to)- 

The function g^j introduced earlier (j3.ip is a sum of four terms. We start 
by separating out these terms. Let 



P+. = — 

(n 



' r=l u Wr s=l Ws s=r+l Ws 



^ d C~l r-1 

^v,j ~ 2Y ^™ "^vwr ^nt»r TT y^^^-^ ~ 

^ ' r=l U liJr j=j + l S = l Ws 



X 

s=r+l uis 

In both expressions above, the first sum is over all values of n 7^ t", and the 
subsequent sums are over all values oi Wr ov Wg u,v. Let further 

^ ' r=l U v' S=l Ws S=T-+1 Ws 

^ ' r=l u v' i=i+l s=l Ws 

X ri E(i-4vi). 

s=r+l tfa 

In the last two expressions above, the first sum is over all values of u 7^ i^, 
the second sum is over all values of v' / u, f , and the subsequent sums are 
over all Ws 7^ u, w'. 
Then 



We define further 'standardised' versions of P^j, Pj, j, Qvj- 
1 

P+ = t Vf Vi-^' 

f ,i _ o^2d Z-,' ™ f"'- uw',. 



(n - 2)2 
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X n E (1 - ^j^i^) n E (1 - ^'^")^ 

s=lws,'w'g s=r+lws,w's 
^ d C-1 

p-. = - Wi^ Vf V Vf , 

X n E (1 - ^^j^o n E (1 - e'l^^)- 

In these two expressions, the first sum is over all values of u ^ v, and the 
remaining sums are over all values oi Wr,w'j.,Ws,w'g ^ u,v. Also, let 



^v,j ^ - 2)2rf E E ^w'^ E E 



X li E (1 - n E (1 - ^v^')^ 

s=lws,'w'g s=r+lws,w'g 
^ d C-1 

^v,j = _ 2)2d E E E ^^y' E E^^"'' 

^ ' r=l u v' i=j+l v" 

X n E (1 - ^:co n E (1 - ^^)- 

In these two expressions, the first sum is over all values of tt 7^ f , and 
the subsequent sums are over all values of v' ,v" ,Ws,Wg 7^ u,v. In these 
standardised versions, the "anchor" nodes u and v for the final products are 
chosen so that the products can be extracted as far as possible as a common 
factor. In the future, we will always use similar conventions regarding the 
ranges of the various summations involved, and the choices of anchor nodes. 
Set 

9v,j = Pv,j + PyJ + Qtj + Qv,j- 

We shall now bound {g^j — g^jl above via upper bounds on the differences 

\P^j Pvjl - P^jl IQtj - QU and IQ- . - Q-J. These differences 
denote the maximum difference of the functions over all load vectors x. 
Noting that, for < j < C and for any n, 

^ V /"i - 1^-?' i^^^' ^ - V f 1 - 1^-?' i-^ ] 

Ws \ J - 

^ E ^Ms - E ^J^i + 
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we see that 

\p:^^-p:[^^\<d\c+if{n-2)<t^\ 

Similarly, 

\P-^-p-^\<d\C + lf{n-2)4>\ 

Also, 

< d{d - 1)[(C + lf{n - 2)(^^ + (C + l)(n - 2)(t>'^] 

^ (n - 2)2<i-i XI XI ^w'' I XI ^S^' {}v'i ~ :^^^zr2 X ^v"v} 

^ ^ r=l u v' v" 

s=lws,w'g s=r+lws,w'g 

d 



< d{d - 1)(C + lf{n - 2)(</.i + </>2) + ^ ^ 

^ ' r=l u 

X n X (1 - i^iJi) n X (1 - ^'t^') 

s=lws,wi. s=r+l WsjW'^ 

^ I ;732 X ^w^'^f'i ~ (n - 2)2 X X ^^"^^ • 

Hence 

< d(d-l)(C + l)2(n-2)(</.V0') + rf(C + l)(n-2)0i 

< 2d^{C + lf{n-2)4>. 

Similarly, IQ"^- - Q-J < 2d^{C + l)3(n - 2)0. 
It follows that 



\gv,j-9v,j\<Gd'iC + l)%n-2)4>, (6.1) | gghat 

and so 

\E[g^j{Xt)] - E[g^jiXt)]\ < U\C + lf{n - 2)E[^(Xt)]. 
In the special case d = 1, the estimates are easier, and we find that 



\gvd - gvA < 4(C + l)(n - 2)0. (6.2) | gghat 1 

Our next task is to estimate |E[^^j(Xt)] — (n — l)gj{zl)\, where (n — 
l)* (^t') is given by 



. E[/,,<,] X] ((n - 1) - E[/,,<,]^ 



2\r-l 



x{{n-l)-nfv,<j-if) 
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C-1 



+ 



(n 



' i=j+l r=l 



x((n-l)-E[A,<,_i]2)' \ 

Here, and throughout what fohows, we abuse notation by writing e.g. E[/^^c'] 
instead of E[/t,.c'(X()]: for the remainder of this proof, ah of our functions 
win be evaluated at Xt. 

We start by estimating the difference between E P^- and 

_ iW+1 2^ (('^ - 1) - E[^,<i] ) ((n - 1) - E[^,<i-l] ) 

r=l 

I ¥j[gvj{Xt)] — (n — is the sum of this and three similar terms. 

Note that 



v,<j,j- 



By (jS.lSp . as n > no(A, d, C, to)i the function 

^ d r-l 



r=l tOr 



= 1 Wa,w'^ 



satisfies, for each t < Iq, 

1 



\fn,v,<j,j{Xt) - E[fu,v,<j,j{Xt)]\ > -y/nlognj < 
By Lemma [^31 as n > no(A, d, C, to)) 



-7 log n 



1 



EKi;/«,-,<i,i] - n^uv] nfu,v,<j,j] < ^V^logn + 8ne-^i°s " < V^logn, 
for each u ^ v, and so 



< \/nlogn < 2^/n\ogn. (6.3) 

n — 2 



Now let Et be the event that \fyj{Xt) — 'Ef^j{Xt)\ < ^^/nlogn and 
\fv,<j{Xt) - E/„,<j(Xt)| < i^^logn for all j G {0,...,C} and v G V^. 

By (|5J5]1 and (|5J6]l . P(;E^) < 16(C + l)ne-^'°s'". 
Note that, on Et^ 



1 



n - 2 
1 



n - 2 
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for each j. Thus, recalhng that 



r— 1 



X 

r=l 



((n - 2f - - l£^)(/.,<, - l^i) 

r=l 

X ((n - 2)2 - - _ 4^-1)) 

we see that, on the difference between fu,v,<jj and 

1 r— 1 

r IE[/.,<,] E[/.,,] J] ((n - 2)2 - E[A,<,] E[/„,<,]' ' 



(n - 2)2rf 



X ( (n - 2)2 - E[/„,<,-_i] E[/,,<,-_i] 



is at most ^d\/nlogn in absolute value. 

Thus, with probabihty at least 1 - 32(C + l)ne-^^°>^ fu. .'u,<j,j is within 
distance ^^/nlogn of E[/„_^^<jj] and within distance ^d?y/nlogn of 



1 r— 1 

X ((n - 2)2 - E[/„,<,-_i]E[/,,<,-_i]) 



It follows that the difference between ^[fu,v,<j,j] and Hu,v,j{t) is at most 
2(f' y/n\ogn. Since, for each u^v and j, | E[/„^<j] — E[/t,^<j]| < (n — 2)(C + 
1) E[(/)2], the difference between '^[fu,v,<j,j\ and Hu,v,j{t) is at most 2d? y/n\ogn+ 
d'^{n - 2)(C + 1) E[(/)2], for each u, v and j. 

Combining the above with (j6.3p . we see that the difference between E[PJ*'^] 
and Hu.v.j{t) is at most Ad'^^/nlogn + d'^n{C + 1) E[(/>2] in absolute value. 

A similar argument shows that the difference between E[P^^] and 

-— -2^E[/,,c]E[/„,,] ^ E[A,i]^((n-2)2-(E[/,,<,])2'''" 



x((n-2)2-(E[/,,<,_i])2) 



d—r 



is at most 4(i2C-y/nlog n + d'^nC{C + l)E[(/)2] in absolute value. 

For E[(5^J, we use an argument identical to the one above, considering 



/ v,u,<j,C - _ 2^2d-l XI XI ™' X v"v 

^ ' r=l v' v" 
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r-l 



s—r+1 Ws,w'g 



to show that the difference between E[Qjj] and 

d 



(n 



nfv,c\ nfv,<j] iE[/.,,] ((^ - 2)' - nfv,<j]) 

' r=l 



X (n-2)^-(E[/,,<,_i])^ 



r-l 



is at most 4(i^y^logn + (fn{C + 1) E[(/)^] in absolute value. 
Similarly, the difference between E[(5~j] and 

C-l d 

■E[/„,c]IE[/.,,] IE[/.,.]E(("-2)'-(K[/-,<*])' 

i=j+l r=l 

xf(n-2)2-(E[/,,<i_i])^ 



(n - 2)2^^ 



is at most A(PC^/nlogn + (PnC{C + l)E[(/)2] in absolute value. 
In summary, we have shown that 

E[g,^,{Xt)]-{n-2)g,iy];) 
< 



r-l 



ngvA^t)] - navA^t)] + ngvA^t)] - (n - 2)g,{yi 

< 8(f{C + lfnE[^{Xt)] + 8d2(C + l)^/nlogn, 

where is the vector with components yt{v,j) = ^^E/„j(Xf). 

Now, each of the components zt{v,j) and yt{v,j) is non-negative, and we 



have Yl^o^tivJ) < 1 and Y!j=Qyt{'"^j) ^ Furthermore, - 

(^^, i) I < for all j. Also, exactly as in the proof of (|7.1|) below, whenever 
y and z are in {x G R'-^+^ : x{j) > for each j, J2j ^U) ^ ^5^}) have 

3 

max - z{j)\. 



n — 1 

n — 2 J o<J<c 
n-l\3 1 



1 



\gk{y)-gk{z)\<3d'{C + l) 
It follows that, for n > 6, 

- 9A^B\ < id\c + If (^j ^73^ ^ 6^'^^ + - 2 ' 

and so, using the fact that < 2d{C + 1) for any load vector x, 

E[5,,,(Xi)]-(n-l)5,(zn 
< S(f{C + l)=^nE[(^(Xi)] + M^{C + l)V^logn + Q(f{C + if + 2d{C + 1). 
As n > (C + 1)2, we may now write 

ngvA^t)] - (n - < 8d'(C + l)3nE[0(Xt)] + 16d2(C + l)V^logn, 
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as claimed. 



□ 



We now need to study the changes of (t){Xt) over time. 

For distinct vertices u and v, and j, A; G {0, . . . , C}, we define 



^u,v,j,k 



n 



k . 

vw' ' 



6^ • 

U,V,J 



n 



;ifu,j 



■ fv,j) 

1 

n — 



n 



V F , 



^ x({n,t!},7i;). 



Then we have (j)'^ = maxu,v,j,k\(l)i^yj^k\^ ^'^ = ^^^u,v,j \(pl^v,j\ and (/>^ = 
maxu,!) (/"i^i,) where all maximisations are over distinct vertices u and v and, 
where appropriate, j,k G {0,... ,C}. These functions are similar to ones 
in [§] : we prove an analogue of Lemma 2 in [3| , leaving some details to an 
appendix, but our task is more complex as we deal with d > 1, and we fill 
in a key point that is dealt with rather brusquely in 0]. 

Once again, our argument uses the discrete chain (Xt)- Let (J^t) denote 
the natural filtration of (Xt). Let Bt-i = {Xg G S for all s < t — 1}. 
For a function / : S ^ M, we define Af{Xt) = f{Xt) - fiXt-i), the 
increment of the function on one step of the discrete chain. Our first goal is 
to provide upper bounds on E[\A^l^^^^f^{Xt)\ \ Pt-i], E[| A02 „ .(Xt)| | Pt-i\ 
and E[| A(/)j^ | in terms of (j){Xt_i)^ valid, on the event Bt_i, for 

all distinct nodes u and and, where appropriate, all j, k G {0, . . . , C}. 

The proof of the following lemma consists of routine but tedious calcula- 
tions, and these are relegated to the appendix. 



lem. appendix Lemma 6.3. Suppose t < Iq and n > no(A, d, C, to)- For p any one of the 



functions 4>\ ^ • ^, (j); 



or 



we have, on B- 



v,j,k> ^u,v,j> '+'u,v' 



t-i, 



ci 



lem . phi-over-time 



where d = 26(1 + l/X)d^{C + l)^ and ca = 64:Xd^{C + l)^. 

If d = 1, we have the same conclusion with ci = 26(1 + 1/A)(C + 1) and 
C2 = 64A(C + 1). 

Now we are in a position to prove the other result required for Lemma [6. 11 
Lemma 6.4. For all t < to, and n > no(A, d, C, to), we have 



E</.(XO < 2e208(A+i)'^^(C^+i)^*o/'^(^^) 



+ 



3 log n 



n 



Ifd 



1, and n > no(A, 1,C, to); we have the improved bound 
E0(Xi) < 2e208(A+i)(C+i)to (^^(^^) + 31ogn^ 
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Proof. Let m = (j){xQ) + and let Ct be the event 

Bt n {(piXs) <m(l + ^y , for all s<t}. 

Let p denote any of the functions 4>\yjk-, 4>1ivj 4't.,v Fo'^ each t, on 
the event Ct-i, we have from Lemma 16.31 that 



E(|p(Xt)|-|/>(Xt_i)| I <E(|Ap(Xt)| I Ft-i)\ <^m[^ + ^) + 

and therefore 



This yields, for each t, 

s=0 ^ ^ 



Cl \* 31ogn C2t 
+ ^ ^ + ^- 



Now let c = 8Ato, and run the discrete chain for en? steps. Note that 

n > no(A, d, C, to) > 10^(A + l)^d^(C + 1)^*^ > (cca)^. 
For t < CTi^, we conclude that 

Cl \ * 2 log n 



E(|p(Xi)|Ic,_0 <m 1 + 



We now show by induction on t that F{Ct) < te"^''''"^^"', for all t < c-n? . 
This is certainly true for i = 0. If the induction hypothesis holds for t — 1, 
then we have 

E(|p(l0lfc) < C^iC^Ti) < Ccn2e-5^'°s^" < ^ 



/n 

and so 



Thus 



logn 



t) <P(5t)+P(a_i) + (C + 2)^n^maxP |p(XO| >E|p(Xt)| + 

where the maximum is over all the functions (puv k j' ^uv j ^^'^ 'Pu,vj noting 
that there are at most (C + 2)^n^ such functions. 
Inequality (|5.5p implies that 

for t < cvP, since the chain starts in 5*1. 
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To establish concentration of measure for the functions | ^ ^ ^ | , I t, j I 
and we proceed as in Sections [H and O Indeed, it is easy to see that 

{n - 2)|K„,,- fc(x)| - fc(y)|| < 2{\\x - y\U + \\x - y\U), 

{n - 2) \(ti,,j{x)\ - < 2(||x - y\U + \\x - y||„), 



and 



(n-2) <liJx)-ctiJy) 



< \\x — y\ 



for all u, V, j and k, and all x and y in S. Calculations exactly as leading 
up to ()5.6p now give, for any of the functions p, any t < an} ^ and any a < n, 



¥{\\p{Xt)\-^\p{Xt 



> 



< 4e 



< 4e' 



-470^ /n 



n-2/ 

Applying this with a = ^^/nlogn gives 



We thus have, using also the induction hypothesis, that 



nct) < \e 



'+(t-l)e- 



as required for the induction step. 

This proves that, with probability at least 1 — have 



ct>{Xt)<m{l + ci/n^f < U{xo) + 



31ogn 



n 



for all t < en?. 

Recall that to = c/8A. Thus, with probability at least 1 — e~^" there 
are no more than cv? events in the continuous-time chain X during the 
interval [0, Iq]. Since 4> is bounded above by C, it follows that, for all t < Iq, 



{xo) + 



31ogn 



n 



< 26^^=1^*0 (^0(xo) + 



31ogn 



n 



Substituting for the value of ci gives the required result, both in the 
general case and the case d = 1. □ 



7. Proof of Theorem 11.11 

sec : proof -main 



We now use the results of the previous section to derive the main theorem. 
We need one routine lemma, showing that the differential equation (jl.Sp is 
Lipschitz with an appropriate constant, in the domain of interest to us. 
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lem. lipschitz Lemma 7.1. Let d and C he positive integers. Let A > 0. The differential 
equation is Lipschitz with constant 8d'^{X + 1)(C + 1)^ on the set {z G 



j^c+i . > Q jgj, each j, z{j) < 1} with respect to the ioo norm. 

For d = 1, the differential equation has Lipschitz constant 2X + 2C + 6 on 
[0, l]*^"*"^, with respect to the i^o norm. 

Proof For < A; < C, 

\Fk{x) - Fk{y)\ < X\x{k - 1) - y{k - 1)| + X\x{k) - y{k)\ 
+k\x{k) - y{k)\ + ik + l)\x{k + 1) - y{k + 1)1 
+\gk-i{x) - gk~i{y)\ + \gk{x) - gk{y)\- 
Now, for X, y G {z G : z{j) > for each j, z{j) < 1}, 

d 

\2x{C)J2x{k)x{< k){l - x{< kfy-Hl -x{<k- if^"' 

r=l 
d 

-2y{C) y{k)y{< m - y{< kfr-\l -y{<k- 

r=l 

<2{2d + d{C + 1) + d{d -1){C + 1)) ^max^ \x{j) - y{j)\ 

<3(i2(C + l) ma^x^|x(i)-2/(i)|. 

Also, 

d C-1 

\x{C)Yx{k) ^ x{i){l-x{<ifY-^{l-x{<i-lff^' 

r=l i=k+l 
d C-1 

-y{c)Y,y{k) ym-y{<ifr\^-y{<^-^?f~^ 

r=l i=k+l 

< 3d^C{C + 1) max \x{j) - y{j)\. 



o<j<c 

It follows that, for A; = 0, . . . , C — 1, 

\gk{x)-gk{y)\<M\C + lfmax\xU)-y{j)\- (7.1) 



So, for < /c < C, for x,y e {z £ R^+'^ : z{j) > for each j, z{j) < 1}, 



o<i<c 

n fnr PF\rh i 



\Fk{x) - Fk{y)\ < {2X + 2C + Qd\C + if) max \x{j)-y{j)\ 



< 8d'{X + 1)(C + 1)^ max |x(i) - y{j)\, 

0<3<C 

and the same bound holds for A; = and k = C. 

For d = 1, it is easy to see that, for /c = 0, . . . , C and x,y £ [0, 1]*^^^, 
\gkix) - gk{y)\ < 6maxj \x{j) - y{j)\, and therefore 

\Fk{x) - Fk{y)\ < (2A + 2C + 6) max \x{j) - y{j)\. 

0<j<C 
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lem . unique-sol 



for each k = 0, . . . ,C. □ 

Lemma 7.2. Let d and C he positive integers. Let A > 0. Let satisfy 
Z]j=o^o(0 = 1 and ^o(^) ^ for all i. The differential equation ( [j.gj) has a 
unique solution (^f) subject to initial condition valid for all times t > 0. 
Furthermore, Ylj=o^i(^) ~ ^ '^"'^ ^t{j) > for all j, for all t > 0. 

Proof. By Lemma l7.ll F is locally Lipschitz with respect to the ^oo norm, 
so the differential equation ()1.3p has a unique maximal solution (^t) valid 
on [0,tmax) for some tmax > 0. Note that Yl'j=o^j(^) ~ ^ ^ ^^'^ ^° 

Yl'j=o^t{j) is constant for all times t < tmaxi and hence is equal to 1. Also, 
F'jiO ^ whenever ^(j) = 0. By standard arguments, S,t{j) > for all j 



and all t < tmax- Thus G {z G 
for all t < tmax, and hence tmax = 



pC+l 



oo. 



> for each \z{j)\ < 1} 

□ 



Lemma 7.3. Let A and to be positive reals, let d and C be positive integers, 
and suppose that n > no{\,d,C,tQ). Let S^o satisfy CoU) ^ for all j = 
0, . . . ,C, and Ylf=o Co(i) = 1- Then, for each v and each t G [0, to]; 

^-nfv,AXt)]-uj) 



sup 

< 



n 



sup 
j 



n 



-jfv,j{xo) - Co(i) 



+ 45AM^(C + l)3U(2;o) + 



3 log n 



n 



,216(A+l)d2(C+l)3to 



For d 



sup 

j 



1, we have 

^—nfvAXt)] 



n 



< 



sup 

j 



n 



-jfv,j{xo) - Co(i) 



+ 45Xto{C + lfi^(j){xo) + 
Proof. For each j, for t < to, we have 



Slogn 



n 



,216(A+l)(C+l)to 



=Co(j)+ / F,iCu)du. 



As before, for u G K and j G {0, . . . ,C}, j) = (n - 1)-^ E[A,j(Xt)], 
and is the vector {zt{v,j) : j G {0, . . . , C}). Also, j4 is the generator of 
the process X. Then, for every j, 

^zt{v,j - 1) - >^zt{v,j) -jzt{v,j) + (j + l)zt{v,j + 1) 
A 



dt 



+ _- -(E[5.,,-i(X0] - IEb.,,(^«)]), 
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so 

zt{v,j) - zo{v,j) 

1 



n — 1 

t 



(^Xzu{v,j - 1) - XzuivJ) - jzuivj) + (j + l)zu{v,j + l)^du 



I Fj{zl)du + \ -n9v,j-i{Xu)]-9j-i{zl)]du 

-^l {:;^n9vAXu)]-gjK))du. 

For each j and t, let et{v,j) = sup„<( \zu{v, j)—^u[j)\, and let be the vector 
with components et{v,j) (j = 0, . . . , C). Let L be the Lipschitz constant of 
the function F, as in Lemma |7. 11 Then we have, using Lemma l6.lt 

||ej||oo < l|eolloo+/ f^lle^lloo + 2Asupmax -E[5^j(X„)] "jdu 
Jo ^ u<t 1 n- I J 

< lleolloo + L I \\el\\^du 
Jo 



By Gronwall's lemma, for each t < t^, 



II < pLto / ||,»^|| 
HWoo ^ I II toll CO 



which gives the required result, since we may take L = 8d'^{X + 1)(C + 1)^ 
by Lemma l7.ll The result for d = 1 follows in an identical manner. □ 

Proof of Theorem 11.11 Let 

^/n,,(xo) - eo(j)| + 45Xtod\C + If [(fixo) + ^^)) 

^ g216(A+l)rf2(C+l)=*to_ 

The previous lemma, along with ()5.19p . yields 

sup|/,,fc(X0-(n-l)6(A;)| > (n - l)u; + log ^ 



w = [ sup 



n 



v,k,t 

< p(sup|/,,fc(Xt)-E/,,fc(Xt)| > V/^logn 
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where the supremum is over ah v G Vn, k £ {0, . . . , C}, and t < Iq. 
The proof of Theorem 11.21 is essentially identical. 



□ 



sec : init-cond 



8. Initial conditions 

Given a vertex v, in order for the functions fv,k{^t) {k = 0, . . . ,C) to 
be well approximated by the solution (^j) to the differential equation (|1.3|) . 
given the initial state Xq of the system, we must choose initial condition 
for (jl.3p in such a way that supj^^^Q ^y eo{v, j) is small, where eo{v,j) = 
l^o(j) — {n — 1)^^ K[f,uj{XQ)]\. For instance, we can take Co{j) = {n — 
l)-^E[f^j{Xo)] for j = 0, ...,C. In addition, there are restrictions on 
allowed initial states Xq, to ensure that (f){XQ) is not too large. 

Clearly, Xq = implies that (j)^{Xo) = 0, so the law of large numbers in 
Theorem 11.11 holds if satisfies Co(0) = 1 and .^o(j) = for j = 1, . . . , C. 

Now consider an initial state obtained as follows. For a constant cq > 0, 
we throw [00(2)] calls onto the network, one at a time. Each call chooses 
endpoints u and v uniformly at random; it is routed onto {u, v} if there is 
spare capacity. Otherwise, it chooses an ordered list of d intermediate nodes 
. . . , Wd) uniformly at random with replacement and is routed onto the 
first route {u, Wi}, {v, Wi} minimising the maximum load of the two links, if 
this route has capacity. If each of the d routes has a full link, then the call 
is lost. Let Xq be an initial state obtained from this [00(2) J -step allocation. 
We will show that, with high probability, (^{Xo) is at most 3n~^^'^logn. 

We start by analysing (j)^{XQ), using the bound 

<p\Xo) = max maxl^- Vli^(Xo)lL(^o) 

u,v:u^v 7,fe \ n — Z 

Y ^ ^ I 



< max max 

u,v:u^v j,k n — 2 



j;ii^(Xo)lL(^o) - EEli^(Xo)lL(^o)] 

to to 



+ max max^-|E[VlL(Xo)lL(Xo)] 

^ to 

- J2 lito(^o) E ^-'(^0) 



+ max max - 

u,v:u^v j,k (n — 2) . 
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Arguments similar to those in Sections|l]and[5]show that li,w{XQ)l'^^{XQ) 

and J2^^u,v^^nw{Xo)J2w'y^u,v'^vw'(^o) are well-concentrated. Specifically, 
there exists a constant 70 > such that for all u, v,j, k, 



n-2 



and also, for all u,j, 



We deduce that, with probability at least 1 - 16(C + ifn'^e-""'^''^^ ''', 
cP\Xo) < max max^|E[^Ii^(Xo)lL(^o)] 



u,v:u^v j^k Tl — 2 

UJ 



n-2 



+ 2 



log n 



n 



(8.1) eq.dev 



From ()8.ip we have that, for sufficiently large n, 

E [fuAXo)fu,k{Xo)] - E fuA^o) E fuA^o) 
= E [{fuAXo) -EfujiXo)) {fu,kiXo) - E/„,fc(Xo))] 

< Qv^logn^ + 16e"'^o^°^'''n2 < ^nlog^n. 

Now we use the fact that, for fixed w, all of the variables Xq({u,w}) are 
identically distributed, to obtain that 

E/^j(Xo)/^,fc(Xo) = (n-l)EIi^(Xo)lL(^o)+(n-l)(n-2)EIi^(Xo)lL(^o) 
for any distinct u and v, and therefore 
1 



{n-iy 



■ E U,, {Xo)f^MXo) - E li^ (Xo)lL (^0) 



< 



n — 1 



Note also that E/u,j(Xo) = (n — l)Eliu)(^o); for any u and tt;. Hence 



E 



lL(Xo)C(^o) - ElL(Xo) eC(Xo) 



< 



1 



+ 



1 1 



n-1 (n-l)2 2 
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log^ n 
n 
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for any distinct u, v, w, and any j and k. The same argument gives 



E 



IL (^o)IlL' (^o) - E (Xo) E IL' (^o: 



< 



log^ n 
n 



i.3) eq.indep-4 



whenever u, v, w and w' are all distinct, and for any j and k. 

Now, by (j8.2p and (j8.3p . for any distinct and v, and any j and k, 



n - 2 



< ^^ + ^1 E[F_(Xo)]E[iLTO] 
n n — 2 I ^—^ 

— EK,(Xo)] IE[lL'(^o)] 



n - 2 



Since all the luw are identically distributed, all the this last term 

is zero, and we have, for any distinct u, v, and any j. A;, 

^|e[^ f,^(Xo)iL™]--^e[^ F,^(Xo) tAx,)] 



n-2 



< 2 



log^ n 
n 



It follows that, for n large enough, 

, log n ^ 



i(Xo) > 3^^) < 16(C + l)2n2e-^«'°s ^. 



Furthermore, 



max max 1 fuj {Xq ) - f^j {Xq ) 

u,v:u^v J n — Z 



< max max 

u,v:Uy^v j n — 2 



\fuAXo)-E[fuj{Xo) 



+ |Aj(Xo)-E[/,,,(Xo)]| , 



and hence, by (jS.ip . 



log n 



(Xo) > 2^^) < 8(C + l)n^e 



2„-7olog n 



For 0^, standard Poisson tail bounds yield that, for each fixed pair {n, v}, 
the probability that there are more than cq log^ n calls with endpoints u and 
V is at most (e/log^ n)*^"'"^ < e~^°'°^ " for sufficiently large n. Thus 

colog^n^ 



{<PHXo) 



> 



n 



43 



Malwina J. LUCZAK 



Hence, as claimed, for n large enough, 

n<p{Xo) > 3^^) < 25(C + l)2n2e-^»i°s'" < e-5Toi°g'". 

9. Extensions 

Theorems 1 1 . 1 1 and 1 1 . 2 1 imply a 'global' law of large numbers approximation 
for the network, that is the number fk{Xt) of links with load k is well 
approximated by the differential equation (jl.Sp . Indeed, for instance, by 
Theorem ll.il when d >2, summing over all the vertices gives the following. 
Let An be the event that, for each k and each t S [0,to]) 

\fk{Xt) - l^yt{k)\ < (sup|/,(Xo) - Qeo(j) 

+ 23(A + l)(to + l)d\C + l)3(nV(^o) + Sn^/^ logn))e2i6(A+i)d2(c+i)3to, 



Then P(^„) < e-2Ti°s In the case d = 1, an analogous result can be 
deduced from Theorem II. 2i It would appear that these results are unlikely 
to be close to best possible: we would expect to be able to approximate 
fk{Xt) with error of order 0{n), up to a logarithmic term, but have not 
been able to prove such a result using our methods. There are several places 
where our argument would need to be improved, including the concentration 
of measure arguments used in the proofs of Lemma 16.41 and of Lemma 16.21 
Our techniques can be adapted to analyse all the other variants of the 
model mentioned in the introduction. More generally, one would expect to 
be able to handle models involving a large system (of size n), where any pair 
of elements (e.g. links) interact at a rate tending to as n — )• oo. (In the 
present model, any pair of links share an arrival stream at a rate of order 
0(l/n).) These extensions may require a modified definition of function (p. 
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Appendix A. Proof of Lemma 16.31 
Proof of Lemma 16. 3[ We start with 4>\^ jj., and note that 



3,k 



n 



2)2 



A E E 



We therefore need to compute and compare the conditional expectations of 
these two increments. 

On the event Bt-\ = {Xg G S" for all s < t — 1}, we have, for any u, v, j, k, 



F 



JKw + U + ^Wut} + H^iw^ - ^iw + 9u,w,j-i - gu,w,j) 



+ iPu,v,w,j-l,k-l + Pu,v,w,j,k)^{Xt-l] 



(A.l) 



+ li_i({^x,t;},«;)(F„^lS^ + Ii+il^+i)(l, 

where, for instance, X{^vwgu,w,j) (-^t-i) is the contribution for the case where 
a call is indirectly routed via the link uw which has load j in Xt-i, and 
v,w,j—l,k—l (Xt-i) is the contribution for arrivals onto the route consist- 
ing of the links {u,w} and {v,w}, with loads j — 1 and k — 1 respectively 
in Xt-i. The term Xt_i({n, u}, t(;)l{ii^][^+^(Xt_i) represents depar- 

tures of calls from the route consisting of links {n, w} and {v, w}, with loads 
j+1 and k+1 respectively, while the term X]u,^„,„ Xt-i{{u, v},w)liwIyy,iXt-i) 
represents departures of calls from the route consisting of links {u, w} and 
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{v,w}, with loads j and k respectively, which, along with the arrivals to 
such routes, have otherwise been overcounted. 

Explicitly, in the expression above, for each u, w, and j, 

q.. ..,, = P+ ■ + ■ + Q~^ ■ + Q~ ■, 

with 

^ d r—1 d 

^u,w,j = C„ _ 2)d ^"'^ X] ^v'uK'L ~ K'u,wJ n ■'^i'i.wj! 

^ ' r=l v' ,Wr s=l s=r+l 

^ d C-1 r-1 d 

^y;w,j ^ (n _ 2¥ ^ ^""^ ^ ^ ^^'"^ H'^"'- ~ 11 ~ 

^ r=l v' ,vfr i=j+^ s=l s=r+l 



^ d r—1 d 

= ^ ^ 11(1 - n (1 - Ji.-'- 

r=l v' ,vfr 5=1 s=r+l 



r-1 

<i-i 

U!,tOs'" 



^ d C-1 r-1 d 

^u,w,j = _ 2)d '^^uw XI ^^w' XI ^^-'u ~ 11 ~ Vty,lj' 

^ r=l v' ,vfr *=i+l *=1 s=r+l 

where, according to our convention, denotes the sum over all v' 7^ u, w; 
in the first two expressions, J2yvr denotes the sum over all wi, . . . ,Wr-i, 
Wr+i, . . . ,Wd such that Ws v' ,u for any s, whereas, in the final two expres- 
sions, denotes the sum over all wi, . . . , Wr-i,Wr+i, . . . ,Wd such that 
Ws ^ v',w for any s. Also, explicitly. 



^ d r—1 d 

Pu,v,w,j,k — / _ y ^ ^uvHiw\w y ] XT^l ~ ^uv,Ws) XI (1 ~ ■'^uifuis ) 



where here denotes the sum over all , . . . , Wr-i, Wr+i, . . . ,Wd such 

that Ws ^ u,v for any s. 

Similarly, on the event Bt-i = {Xs G S for all s < i — 1}, 



= {( E ^-') E [-fliy^ + u + mi'+m-^'-K.. 

+ gu ,w,j—l 9u,w,j 



+ ^(^W ~ ^L' + 9v,w',k-l - 9v,w'>k)\ 
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+ lt_i({n,t;},«;)(Ii^lL + liril'i')(^t-i)- (A.2) 

Comparing corresponding terms in the two expressions (jA.ip and (IA.2P 
gives that, on the event Bt-i, 



< {2j + 2k + 2 + 4:X)(j)\Xt^ 



~^ ^7 0\2 {Pu,v,w,j-l,k-l + Pu,v,w,j,k){^t-l) 

[n — Z) 

~l~ \(^u,v,j,k—l ^u,v,j,k ~l~ (^v,u,k,j—l ^v,u,k,j\{-Xt—l) 1 (-A'^) 

where 



n-2 

w 



n-2 

Bounding the middle two terms in the expression (|A.3p above is straight- 
forward: note that each Pu,v,w,j,k is at most d/{n — 2), while 

Xt.i{{u,v},w)<{n-2)(t>^{Xt-i), 

and so, on Bt-i, 

a(2) + L6a(2)j)e[|A<,,,,,(X,)||J-,„i] 

< (4C7 + 2 + 4A)(/>i(lt_i) + ^ + 203 

n — 2 

~l~ |Oti,i;,j,fe— 1 0,u,v,j,k ~l~ ^v,u,k,j — l Q-Vjnjfcjj I (^t— 1 ) 

< (4C7 + 4 + 4A)(/)(Xt_i) + ^^ 

n — 2 

~l" \^u,v,j,k—l 0,u,v,j,k ~l~ Oi^^u,k,j — l 0,v,u,k,j\{Xf^i'j . 

We will now show that both 
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are close to their 'standardised' version 

d 



u,v,j,k 



2)2d+2 X] X] ^w"' X] ^^'f' ^2 X] 

' r=l w w' v' v" 



X n E (1 - ^iMo n E (1 - ^i^'^') 

s=lws,w'g s=r+l WsyW'g 

+ u.r TD- r>+ 



Analogous bounds hold if i^J^^j is replaced by Pu,w,j^ Qu,w,j^ °^ Qu,i 
First, an elementary calculation similar to earlier ones shows that 



— V P+ t -P+ 

_ 9 u,w,j vw u,v,j,k 



n-2 , 
< d{d-l){C + lf4>^ + 

r-l 



1 



(n-2) 



<i 
■it) 



r=l w 



n E (1 - ^1^.1^^) n E (1 - ^v^') - 

s=r+l Ws,w'g 

2 



< d{d-l){C + if[<p' + (ti' + 

d 
r=l 



n-2 
1 



n-2 



n-2 

^ I (n - 2)2 E ^uw^vw E ^^"^'1 

~ (n - 2)4 E E E ^^''■^ E ^^""ly 

w^u,v w'^u,v v'^u,w v"^u,w 

and similarly that 

r)^2 / y uto-^i)«i / ^ v'u v'w 



in - 2y 



1 



_ 2)4 Z^ Z^ V'«)| 

WJ^U,V w'y^U.V V'y^U,W V"y^U,W 

< (^ + l)'^' + |7;rr2F ^ ^""'^''^ ^ ^ 

~ (n — 2)4 E ^"w' E ^^w)' E E ^v'^w\ 



w^u,v w'^u,v v'^u,w v"j^u,w 



(C + !)(</>! 
< (C + l)(2^^ + ^^ + ^) 



n-2 n-2 



n-2 



) 



n-2. 
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It follows that 



n-2 

w 



+ d(C + l)(2</,i+<^2^-^) 



n-2 
2 

n-2 



If d = 1, we may replace the above bound by (C + 1)(2^^ + + ;^)- 
Similarly, but slightly more easily, 

Hence 

Similar calculations for Pu,w,ji Qu.w,j^ Qu,w,j show that, for each u, v,j, k, 

K,v,j,k\ < 4Ad2(C + 1)3 (2cj>^ + <^2 ^ . 

If d = 1, we have the improved bound \au,v,j,k\ < 4A(C + l){2(f)^ + (f)'^ + ;^). 
Hence, on Bt-i, we have, for all u, v, k,j, 

(K2) + L6^(2)0^[|^<-^>^^*)|'-^*-^] 

< (4C + 4 + 4A)</.(lt_i) + ^ + 16Xd\C + If (3ct^{Xt-i) + 

n — 2 \ n — 2/ 

< 52(A + l)d'{C + ifHXt-i) + ^^^^'(^ + . 
For d = 1, we have the improved bound 

46A(C + 1) 



< 52(A + l)(C + l)</.(Xt_i) + 



n-2 



Next, we consider ^(pl^^j = ;^A(/„j - f^j), for G y„ and < j < 
C. We have 

< {Mfu,j-i - fv,j-i\ + (A + - fv,j\ + (i + - fv,j+i\ 

+ M9u,j-i — 9v,j-i\ + M9u,j — 9v,j\^iXt-i). 
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We now refer to the functions Puj^ uj^ Quj ^^'^ Quj defined in the 
proof of Lemma 16.21 as well as their sum Quj. Using (j6.ip . we have that 

aQ + [6x(^^\)n\HfuAXt)-f.AXt))\ I -^i-i] 

< ^{2X + 2C + l){n-2)(j)^ + 2A(f{C + lf{n-2)(l) 

+ X\guj-i — gv,j~i\ + Mduj — gv,j\'^i^t-i)- 

Ud = 1, the term 24(i2(C+ - 2)0 becomes 16(C + l)(n- 2)(/), by 1^ . 
An easy calculation shows that, for n > 4, and each u, v and j, 

I U,J V,]\ 

^ d r-1 

(n - 2)2f^ ^ ^ ^ ^""'r^M'l); n X] ~ 

X n E (1 - - 7^7^ E E ^u'. E 

s=r+luis,ui3 r=l u' Wr,w'^ 

X n E (1 - ^LMo n E (1 - ^^^^") 

s=lws,w'g s=r+lws,'wi. 

< d^c({n - 2)^2 + 2) + X: E ^u'u E 



in 



r-l 



n E(i-it^.iia) n E(i-ii^'Ji^r') 



s=r+l ■ 



r=l u' Wr,w'^ 



r-l 



X 



n E (1 - ^^Mo n E (1 - lic'ii::^'^ 



< d\C + 1) (^(n - 2)<^' + 2j + (i(2(i + 1)(C + 1) ( (n - 2)(/)^ + 2 

< 4d2(C + l)(^(n-2)02 + 2). 
Similarly, for each n, f and j, 

|P- ■ - 1 < U\C + 1)2 ((n - 2)</)2 + 2 ) ; 



IQ:,, - QX^,\ < M'{C + 1) ( (n - 2)</)^ + 2 ) ; 
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\Q-j - Q-j\ < 4d\C + If ((n - 2)<t>' + 2) . 
It follows that for n> A, and each u, v and j, on Bt-i, 



+ [6A 



n\Aci>i,j{Xt)\ \ ft-i] 

< (2A + 2C7 + 1 + 24d2(C + l)^ + 32\(f{C + lf)(l}{Xt-i] 
64 o , _ 

+ 

n-2 

3 64:X(f{C + lf 



n-2 ^ 

< 35(A + l)d2(C + l)30(Xt_i)+ _ 
In the case d = 1, we may obtain 



n 



+ [6A 



< 35(A + l)(C + l)(/.(lt_i) + 



E[|A<,,,.(X,)||7;_i] 
64A(C + 1) 



n-2 



Finally, we consider the expectation of the absolute value of A(^^^^(Xt), 
conditional on Tt-\- We have 



+ [6A 



E 



( 



AQ) + L6A(^)J 



n-2 



)e[a( ^ 



~ E *-.({«,*»)+;S(i-^En{i-c»-' 

w^u.v ^ ' w s=l 



on event -Bt_i. So we have 



n 



E 



t-i 



< 



n-2 



Xt-ii{u, v}, w) + Xlg < 4>\Xt-i) + 



n-2 



For p any of the functions under consideration, we now have 



n 



+ [Q\ 



\Ap{Xt)\ I Pt-i] 



< 52(A + l)d2(C + lf.^(lt_i)+ _ 

For n > no, A(2) + L6A(2)J > 6A(2) > 2n^ , and the lemma follows. 
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